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The energy-momentum and angular momentum contained in a spacelike two-surface of spherical 
topology are estimated by joining the two-surface to null infinity via an approximate no-incoming- 
radiation condition. The result is a set of gauge-invariant formulas for energy-momentum and 
angular momentum which should be applicable to much numerical work; it also gives estimates of 
the finite-size effects. 
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I. INTRODUCTION 

The exciting successes in recent years of numeri- 
cal treatments of general-relativistic systems have given 
point to longstanding theoretical challenges involved in 
explicating Einstein's theory. Numerical relativists are 
now able to calculate with good accuracy solutions to 
the field equations. But how is one to extract from all 
these data the key physically significant quantities? 

Energy-momentum and angular momentum would be 
among the most important, and the strongest hope is 
that good definitions of these will exist at the quasilocal 
level, that is, on generic acausal two-surfaces (perhaps 
restricted to have spherical topology) . If such definitions 
do exist, and if, moreover, they can be formulated in 
such a way that it is possible to meaningfully compare 
the quasilocal quantities associated with different sur- 
faces, they would be a powerful tool for understanding 
the exchange of energy-momentum and angular momen- 
tum between strongly generally relativistic astrophysical 
systems. At present, however, the problems involved in 
developing a quasilocal kinematics seem so profound that 
one is driven to look for more modest approaches which 
will still provide useful results in broad categories of cases 
of current interest. 

The goal of this paper is to get working approximations 
for the total energy-momentum and angular momentum 
of isolated systems which are suitable for contemporary 
numerical use. Most numerical codes only evolve the sys- 
tem throughout a finite volume of space-time, so one 
has available data on large but finite two-surface S. I 
shall show that there is a reasonable "poor man's" no- 
incoming-radiation condition which can be used to ex- 
trapolate the data on § to future null infinity 3 + , where 
the Bondi-Sachs definition of energy-momentum and 
twistorial treatment of angular momentum apply. It 
turns out — nontrivially — that this approach allows 
for a comparison of the energy-momentum and angular 
momentum as S is moved forward in time; one thus has 
measures of the energy-momentum and angular momen- 
tum emitted in gravitational radiation. 

While a certain amount of work is required to derive 



these measures, the result is a compact set of gauge- 
independent formulas which should be usable by numer- 
ical relativists. 



A. A No-Incoming-Radiation Condition 

The state of a general-relativistic system is specified 
not only by its material degrees of freedom but by its 
gravitational data. In principle, for most modeling one 
would like to fix those data to correspond to no incoming 
radiation. 

To implement the no-incoming-radiation condition ex- 
actly would be very difficult. (One would have to solve 
a hard inverse problem, finding what constraints on 
Cauchy data led to the required behavior of the solution 
in the distant past. Numerical workers typically spec- 
ify the data as well as they can before there are strong 
interactions, and then discard any early transients as po- 
tentially due to spurious incoming radiation.) However, 
this is a practical difficulty and does not affect the va- 
lidity of the condition as the correct restriction on the 
data. 

We may make use of this observation at the two-surface 
§ (assumed to be large and approximately spherical), as 
follows. Consider the null hypersurface Np^ys orthogonal 
outwards from §. Radiation incoming to the future of 
§ would leave its profile on this surface, the transverse 
component being measured by the Weyl component 
in the standard Newman-Penrose formalism. (One has 
= ^ abcdo a o b o c o d , where o A is a tangent spinor 
to Nphys, that is, one has l AA = o A o A for l a the null 
tangent to 3Sf p h ys . ) We may thus take as an approximate, 
"poor man's," no-incoming-radiation condition that 
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$o = on the null hypersurface outgoing from § (1) 

(and there is no matter crossing this hypersurface) . More 
precisely, we consider embedding § in a space-time with 
the same first and second fundamental forms, but we 
discard N p h ys and replace it with a null hypersurface INT 
(orthogonally outwards from S) with "to identically zero 
on 3Sf. (Note that condition (Q]) is imposed only on the 
single hypersurface Ji determined by §. We shall take up 
below the question of what to do when § is evolved.) 

The spin-coefficient equations that effect transport 
outwards from § along Ji can be integrated, so one can 
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work out the asymptotic behavior of the fields as one ap- 
proaches null infinity on §. One can then use the Bondi- 
Sachs energy-momentum and the twistorial angular mo- 
mentum. 

A number of remarks should be made about this. 

First, according to the Sachs peeling property, it is the 
component has the most rapid fall-off along outgoing 
null geodesies (^o ~ r ~ 5 where r is an affine parameter), 
and so one has good reason to think that if § can be 
regarded as in the asymptotic regime, then 5*0 = should 
be a good approximation. 

Second, another way of viewing the nature of the ap- 
proximation is that, because Nphys is at a finite location, 
the incoming waves need not be exactly transverse. That 
is, the extent to which ^0 = fails to implement the 
exact no-incoming-radiation condition is the extent to 
which incoming waves would not be transverse at 3Nf p h ys - 
(This means in particular that the approximation might 
be a poor one if S were substantially wrinkled. How- 
ever, in numerical work one typically has surfaces § which 
are very nearly spherical. If one did have to deal with 
substantially wrinkly surfaces, one could generalize the 
present method by considering an outgoing null congru- 
ence other than the one orthogonal to S, adapted to the 
ambient geometry.) 

Finally, one might consider using the asymptotics of 
the field along N for waveform extraction; cf. e.g. 
(Since these asymptotics only give the field at an in- 
stant of retarded time determined by S, to really ex- 
tract a wave-form one would need to apply the proce- 
dures here for an evolving family of surfaces §(r/) with 
associated null hypersurfaces N(r]), for r\ in some in- 
terval J.) Such a procedure would not be exact, of 
course; its validity would be limited by the applicabil- 
ity of (TT]). While this is certainly natural, it raises is- 
sues beyond those treated in this paper. This is because 
the points which are presently problematic in the case 
of wave-form extraction and in the case of computation 
of energy-momentum and angular momentum are differ- 
ent. For energy-momentum and angular momentum, re- 
solving gauge ambiguities is essential, and condition ([T]) 
allows us to do this (as a mathematical procedure, ir- 
respective of the degree to which the condition accu- 
rately models the physical space-time); for wave- forms, 
the gauge choices for the observer are trivial. Thus the 
definitions here allow one both to compute the energy- 
momentum and angular momentum, and to estimate the 
finite-size effects involved, from data on surfaces §(77), 
but the problem of wave-form extraction necessarily in- 
volves additionally the question of how accurately condi- 
tion models the geometry of the physical space-time. 
One should also bear in mind that, for present work, 
computations of (say) the energy-momentum or angular 
momentum emitted in gravitational waves which were 
accurate to a few per cent would typically be quite ad- 
equate; whereas extraction of physical information from 
wave-forms may require rather more accurate modeling. 
Because of this sensitivity, if the ideas used for asymp- 



totics here are to be helpful to the wave-form extraction 
problem, they will probably be so when combined with 
other physical insights. One should probably use data, 
not just on the surfaces §(rj), but from the portion of the 
physical space-time interior to them, to extrapolate the 
wave-forms. 



B. The Role of Idealization 

The procedure used here turns on embedding the two- 
surface S of interest in a mathematically constructed 
space-time and then evaluating that space-time's kine- 
matics. One might think at first that this is less desirable 
than actually working out the asymptotics along 3Nf p hy S 
to infer the Bondi-Sachs energy-momentum and twisto- 
rial angular momentum in the physical space-time. For 
most purposes, this is not the case, however. 

The point is that the Bondi-Sachs treatment applies 
to idealized isolated systems; for practical purposes one 
must choose which portion of a real (or numerical) sys- 
tem is to be regarded as the isolated component, and 
there may be several such choices. Consider, for in- 
stance, a system which at several widely separated in- 
tervals emits bursts of gravitational radiation. Each of 
these bursts contains, not just the outgoing transverse 
wave front, but smaller trailing, non-transverse, pieces; 
also each of these bursts will generate a certain amount 
of back-scatter via nonlinearities. If one wanted to truly 
be in the mathematically exact Bondi-Sachs asymptotic 
regime for this system, one would have to go outward so 
far along N p hy S that one had passed any slight trailing 
fields and backscatters due to very early emissions — in 
principle, one would have to know the history of the sys- 
tem in the arbitrarily distant past to do this. One rarely 
wants to do this; one would rather think of the system 
as to good approximation isolated in the interval around 
one burst — and if necessary then worry about the fact 
that the isolation is not perfect. 

Thus in most situations the task is not to construct 
the null infinity and kinematics of the entire space-time, 
but to determine how to measure the kinematics of a 
large but finite accessible region. On the other hand, 
the main reason for considering null infinity (together 
with the Bondi-Sachs energy-momentum and the twisto- 
rial angular momentum) — that it provided an invariant 
treatment of quantities of interest — remains valid. So 
our approach is based on constructing a strictly well- 
defined null infinity from the data on the finite surface 
S, in order to have an invariant energy-momentum and 
angular momentum. These are the energy-momentum 
and angular momentum which would be ascribed to S, 
were it in a space-time satisfying (fl}. 

The question of whether one really is in a regime which 
satisfies the Bondi-Sachs asymptotics to a given approx- 
imation — which is the question of whether the energy- 
momentum and angular momentum constructed here are 
stable when S is moved outwards along Nphys — is im- 
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portant but lies beyond the scope of this paper, for one 
cannot investigate this stability based on data at a single 
surface S (nor from a finite family of §(77) 's). 

C. Contemporary Numerical Work 

Numerical relativists who work with codes aimed at 
treating generic space-times have recognized the impor- 
tance of invariant and theoretically justified measures 
of energy-momentum and angular momentum. Indeed, 
this concern plays a part in the choice of some groups 
to use characteristic (or mixed Cauchy-characteristic) 
codes. Such codes, if cast in a form admitting a clean 
extension to null infinity, ought to allow the extraction 
of the Bondi-Sachs energy-momentum and twistorial an- 
gular momentum. 

Characteristic (or mixed) code computations of the 
Bondi energy (often referred to in the literature as the 
Bondi mass, for historical reasons) have been done in a 
number of cases, although typically they are hampered by 
the fact that the codes themselves are not usually cast in 
Bondi coordinates, so that non-trivial gauge transforma- 
tions are required. These also make the computation of 
the Bondi momentum and the angular momentum diffi- 
cult. (See ref. [2] for a review.) The results here may help 
streamline such computations, for the formalism to be 
developed automatically produces the required gauges. 

A great body of numerical work, however, is based on 
"3 + 1" formalisms rather than characteristic or mixed 
ones. For these formalisms, the problem of gauge invari- 
ance for energy-momentum and angular momentum has 
been more severe. 

Most contemporary attempts to extract information 
about the total energy-momentum in the 3+1 formalisms 
can be usefully thought of as based on the Bondi-Sachs 
energy-momentum loss formula [16]. (They have often 
been justified by other means; however, the Bondi-Sachs 
formula would be the broadest and most theoretically 
secure starting-point.) This formula identifies the rate 
of change of the Bondi-Sachs energy-momentum, with 
respect to Bondi retarded time, as an integral of the 
squared modulus of the "news function" with respect to 
the measure induced on an asymptotically large sphere 
by the Bondi coordinates. Thus to use this formula to 
recover the radiated energy-momentum one must know: 
(a) that one is in the asymptotic regime; (b) the Bondi 
coordinate system, both to identify the measure on the 
sphere correctly and to do the integral over retarded time; 
(c) the news function. (The news function can be given 
as the integral of a curvature component — the compo- 
nent depending on the Bondi coordinates — with respect 
to Bondi retarded time.) 

Numerical work in the 3 + 1 formalism has not yet 
implemented any systematic transition to Bondi coordi- 
nates. Thus what is actually done is to use the numer- 
ical angle and time coordinates to compute the curva- 
ture component and integrals required for the energy- 



momentum- loss formula and its integral (e.g. [3J, |4j|). 
Consistency checks are then done by studying the sta- 
bility of the result as the extraction radius is increased. 
However, the lack of gauge control makes it impossible to 
know what these numbers really signify. The stability of 
the actual, gauge-invariant, energy-momentum could be 
either better or worse than the cited numbers, depend- 
ing on whether the gauge freedoms exacerbate or mask 
extraction problems. 

The present approach overcomes the concerns about 
gauge by giving formulas for the Bondi-Sachs energy- 
momentum and its evolution in terms of gauge-invariant 
quantities on the extraction surface §. As far as the ques- 
tion of the stability of the results with increasing extrac- 
tion radius goes, this is an issue which one can only inves- 
tigate directly, by considering larger and larger surfaces. 
However, the present work does give one the confidence 
that in such an investigation other potential error sources 
have been controlled. 

The situation for angular momentum has been more 
difficult than for energy-momentum. In the first place, 
there has been for some time no really theoretically satis- 
factory formula; and in the second, the angular momen- 
tum depends on curvature terms deeper in the asymp- 
totic expansion, which are still more sensitive to the cor- 
rect choice of Bondi frame. The approach given here 
overcomes these difficulties. We use the recent twistorial 
definition of angular momentum, which appears to be 
theoretically satisfactory. We identify the Bondi frame 
exactly on N, and the finite-size contributions may be 
read off directly from the formulas here. 



D. Recent Theoretical Work 

The problem of estimating kinematic quantities in 
terms of numerical data has been taken up by two sets 
of authors recently. 

Gallo, Lehner and Moreschi [5| (see also |6|) raised 
many of the concerns motivating the present work. They 
emphasized the importance of extracting invariant in- 
formation, and also of considering finite-size effects. 
They gave an approach to estimating the Bondi mo- 
mentum which (while presented somewhat differently) 
can be viewed as assuming that the two-surface S is 
only infinitesimally separated from J + , and computing 
the Bondi-Sachs energy-momentum at the correspond- 
ing cut. (While it may seem odd to speak of a two-surface 
infinitesimally separated from infinity, it is a well-defined 
concept from the point of view of the conformally com- 
pleted manifold, and amounts to assuming that one only 
needs leading terms in the appropriate asymptotic expan- 
sions.) Their results therefore correspond to a limiting 
case of some of those here. 

Deadman and Stewart 0] recently discussed the esti- 
mation of the Bondi-Sachs energy from numerical data. 
Their approach is rather different; it is based on con- 
structing a transformation from the coordinates of the 
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numerical evolution to Bondi-like coordinates. As these 
formulas are deduced by considering finitely many terms 
in the asymptotic expansions and extrapolating, this 
work could also be viewed as based on the notion that § 
was infinitesimally separated from 3 + . 



E. Some Technical Points 

A significant feature of the present approach is that, 
with it, all the quantities of interest can be expressed 
in terms of standard Geroch-Held-Penrose (the boost- 
weight-covariant version of the Newman-Penrose calcu- 
lus) quantities at §, and the energy-momentum and an- 
gular momentum are given as natural integrals of these 
at S. The formulas derived here include within them- 
selves all necessary changes to refer to Bondi frames; no 
separate computation of Bondi coordinates is necessary. 

The question of how accurately the Geroch-Held- 
Penrose quantities can be computed at § of course de- 
pends on the particular code. Presumably the most diffi- 
cult one to measure accurately is ^i, which is central in 
computing the angular momentum (and also, because of 
finite-size effects, contributes to the energy-momentum). 
It should be emphasized that the question here is only 
that of the computation of at §, not of its inferred 
asymptotic value (an issue raised by Deadman and Stew- 
art 0). 

A second point is that we need not take up the delicate 
questions of just what degree of smoothness or peeling is 
encoded in the numerical solution. This is because we 
have separated the question of computing the energy- 
momentum and angular momentum from the question of 
finding their limiting values at 3 + : our results are given 
entirely in terms of data at S. 

The approach here also allows one to quantify finite- 
size-effects, and so provides a useful consistency-check 
on the degree to which § is "effectively at infinity" (using 
only data at §). The integrals for the energy-momentum 
and angular momentum at null infinity are given in terms 
of the asymptotic values of the curvature quantities $1, 
^2, ^3; here, those asymptotic values appear as the val- 
ues at § (suitably scaled) plus correction terms. Those 
correction terms, then, are a measure of how removed S 
is from null infinity. We also, as importantly, are able to 
quantify how strongly the inferred structure of null infin- 
ity is subject to finite-size effects as evolution proceeds. 

In numerical work, the surfaces § are typically large co- 
ordinate spheres, and their first and second fundamental 
forms appear as slight perturbations of the values they 
would have for large spheres in Minkowski space. In par- 
ticular, the convergence p§ of the outgoing congruence 
is a slight perturbation of —R^ 1 (by 0(R~ 2 ) or less), 
and the shear cr§ is expected to be 0(R~ 2 ) if there is no 
incoming radiation; here R is the radius of the sphere. 
Thus to good approximation 

|o-s|<|ps|. (2) 



(If there were incoming radiation, one would expect <j$ 
to go like a dimensionless number — the news function 
of that radiation — over R. Thus a small amount of in- 
coming radiation would not upset this inequality.) This 
is helpful, for the asymptotic forms of VP3 on N are 
simplified in this case, and we make this approximation 
in computing them. If more accuracy is needed for par- 
ticular work, probably the most efficient approach would 
be to compute ^2, ^3 perturbatively in crg/ps to the 
required order. (In such computations, note that while 
fs| <C \ps\, the angular derivatives of ps and erg may 
very well be of the same size. Thus one must be careful 
not to discard at one stage terms whose angular deriva- 
tives may be essential later.) Curiously, while ^2, ^3 are 
very complicated, all other elements of the calculation 
are manageable; in particular, even the exact forms of 
the asymptotic spinors and twistors are simple. 



F. Evolution 

I have so far emphasized that, by casting the problem 
of measuring energy-momentum and angular momentum 
at a finite two-surface in a certain form, one can defer 
the difficult questions surrounding quasilocal kinematics. 
However, these questions must be faced to some degree 
when we consider the evolution of the system, and the 
question of how to compare the energy-momenta and an- 
gular momenta at two different two-surfaces. This is be- 
cause the auxiliary space-times constructed via the poor 
man's no-incoming-radiation condition from the two sur- 
faces are not the same, and so it is not obvious how to 
identify the spaces on which their energy-momenta and 
angular momenta take values. Indeed, we must antici- 
pate on physical grounds that unless the extraction sur- 
faces are large enough there will be no way of identifying 
their auxiliary null infinities which preserves all of the 
usual structures. 

This is an instance of a more general problem for 
quasilocal kinematics: how is one to compare the kine- 
matic quantities associated with different two-surfaces? 
Quasilocal kinematic proposals are not well-enough de- 
veloped at present to take up this problem, but the de- 
gree to which quasilocal kinematics will be useful depends 
very largely on the degree to which it can be solved. 

It turns out that, for us, there is a natural approach 
to this problem which fits well with structures previously 
developed for treating angular momentum at null infin- 
ity. We shall see that there is a natural way to compare 
the null infinities from two surfaces §, §' infinitesimally 
separated in time; this procedure can then be integrated. 
Potential finite-size effects show up in that the identifi- 
cations of the null infinities are not via Bondi-Metzner- 
Sachs transformations, but via more complicated mo- 
tions, unless the extraction surfaces are large enough. 
Thus one can say that while a single null infinity is being 
used, to the extent that finite-size effects are important 
the null infinity has a weaker structure than is conven- 
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tional. 

More precisely, we suppose we have a one-parameter 
family of two-surfaces §(77) (for 77 in some interval J) 
foliating a timelike three-surface T, with 77 increasing to- 
wards the future. Each two-surface §(77) is embedded 
(with the same first and second fundamental forms) in 
an auxiliary space-time M(rf) defined by taking \I/o = 
along a null surface N{r)) orthogonally outwards from 
§(77). Then the constructions already described give a 
null infinity J + (77) for each M{rf), with the null geodesies 
orthogonally outwards from 8(77) defining a preferred cut 
6(77) C 3 + (t7). Thus really we have a bundle of space- 
times {M(rf) I 77 e J}-[1Z| (We do not have a single 
space-time for which the condition ^0 = holds on a 
local foliation of null surfaces; that would generally be 
far too restrictive a condition to impose.) 

With these structures, it turns out that there is a nat- 
ural way to identify the null infinities J + (77) for the differ- 
ent 77's. The key step is to make the identifications at an 
infinitesimal level; one can then integrate. The main issue 
then comes down to understanding how one should define 
the cut of 3 + (j]) corresponding to the null geodesies or- 
thogonally outwards from the two-surface §(rj + drj), that 
is, for ri + di] infinitesimally differing from 77. Once this is 
done, one can fix both the identification of the generators 
(because one gets a point-to-point mapping of the cuts) 
and the supertranslation freedom. 

The actual identification we need arises from natural 
isomorphisms. Let us begin with a vector field w a in 
-Wphys connecting §(77) to §(77 + drj). (That is, the field 
w a is tangent to T and w a V a ri = 1.) At each point of 
§(77), we may consider the Jacobi field along the outgo- 
ing null geodesic in M p h ys whose initial value is w a at 
that point (and whose initial velocity is chosen to make 
the field represent a null geodesic outwards orthogonal 
to §(77 + drj)). We thus get a family of Jacobi fields, 
over §(77), which represent null geodesies orthogonally 
outwards from 8(77 + ^77). All of this so far is in M p h ys - 

Now, any Jacobi field is specified by its initial data. By 
the construction of M(rj), there is a natural isomorphism 

of the tangent bundles T(M phys ) = T(M(r})) 

(See footnote [17]. The symbol T which occurs as part of 
the notation T(X) the tangent bundle of X should not 
be confused with the isolated T representing the timelikc 
three-surface foliated by the S(?7)'s.) Thus we may natu- 
rally identify the Jacobi fields we found above with fields 
in the auxiliary space-time M(rj). The limiting values of 
these fields at 'J + (n) we take to define the displacement 
of the cut corresponding to §(77 + drj) from 6(77). (This 
will be well-defined independent of questions about the 
asymptotic behavior of the original Jacobi fields in the 
physical space-time.) It is this definition we required 
for the identification of 3 + (f]) with J + (77 + drj), and its 
subsequent integration to give identifications of the null 
infinities J + (t7') for different values of 77'. 

As noted above, these identifications will not be per- 
fect Bondi-Metzner-Sachs motions, because of finite-size 



effects. The formulas we derive for the evolution of 
energy-momentum and angular momentum apply even 
in this case. However, for purposes of extracting the to- 
tal energy-momentum and angular momentum, substan- 
tial finite-size effects (that is, non-Bondi-Metzner-Sachs 
identifications) should be regarded as signaling that the 
extraction surfaces have not been taken to be distant 
enough. The finite-size results are more of interest at 
present in that they may provide clues about how to de- 
velop quasilocal kinematics generally. 

Having discovered that there is a well-defined null in- 
finity (if with somewhat weaker than usual properties) 
for the family of extraction surfaces {§(77) | 77 6 J}, it is 
natural to ask if one cannot construct a single asymptotic 
regime for the physical space-time M p h ys to which this 
null infinity is attached? In some sense, this is provided 
by the bundle {^(77) 77 G J}, which can be attached 
to -Mph ys along T. However, this bundle is not usually 
a space-time (it does not admit a metric structure com- 
patible with the geometry of the N(?7)'s, since the con- 
dition ^0 = is generally too strong to impose on a 
foliating family of hypersurfaces). While such construc- 
tions might be of some interest in the general problem of 
defining asymptotic regimes, in this paper there will be 
no reason to make use of them; considering the individ- 
ual auxiliary space-times M(r}) and the identifications of 
their null infinities will be what is relevant. 



G. Outline and Conventions 

The next section of this paper outlines the integra- 
tion of the Newman-Penrose equations under the "poor 
man's" no-incoming-radiation condition. In Section III, 
the asymptotic reference frames of Bondi and Sachs are 
introduced. Section IV gives the computation of asymp- 
totically constant spinors and the kinematic twistor. Sec- 
tion V gives the Bondi-Sachs energy-momentum, and 
Section VI the twistorial angular momentum. Sec- 
tion VII derives the formulas for comparing the energy- 
momentum and angular momentum (as well as the defor- 
mation of the numerical coordinates relative to the Bondi 
coordinates) as the two-surface is moved forward in time. 
Section [Villi is a Users' Guide to the results; it summa- 
rizes what procedures and equations would be needed in 
numerical work. 

A reader wishing simply to use the results here can 
use Section IVIIII as an index to the paper. (Since the 
twistorial angular momentum is new, the reader wishing 
to use this will probably want to read the introduction 
to Section VI and also VI A.) 

All material necessary for understanding this paper 
and not otherwise cited can be found in ref . [8j , whose no- 
tation and conventions are used. Ref. Q gives the twisto- 
rial treatment of angular momentum at J + ; ref. [Io| is an 
account of it for non-specialists. 

The standard literature uses the symbol A for three 
different things: a spin-coefficient, a rescaling factor for 
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the spinor dyad, and the angular potential for the shear. 
We shall use Ab for the angular potential for the shear 
and A for the rescaling o A H> Xo A . (The spin-coefficient 
will be — a'; it will play little explicit role in this paper.) 

The metric signature is H — . The symbol ~ will 

be used to denote asymptotic equality as one moves out- 
wards along null geodesies. The symbol « stands for 
equality modulo o( | erg 1/ |ps | ) • All logarithms are natural. 



II. ASYMPTOTICS FROM LOCAL DATA 

The aim here is to compute those asymptotic quanti- 
ties we will need — enough to find the energy-momentum 
and angular momentum at J + along N — in terms of lo- 
cal data on §. This can be reduced to a series of integra- 
tion problems. Some of the details of the (lengthy, but 
straightforward) integrations are omitted. 

What we require are the asymptotic forms of the cur- 
vature quantities ^i, ^2, ^3, as well as the operators 8, 
3' and the shear a. (We also must verify that certain 
other spin coefficients have an appropriate asymptotic 
decay, even if we do not use their values; however, these 
decays will be obvious from the general forms of the in- 
tegrals determining them, except in one case, which we 
shall treat explicitly.) And we must identify the correct 
Bondi frame, that is, we must be sure that when we take 
the asymptotic limits, by looking at two-surfaces reced- 
ing to infinity in null directions orthogonally outwards 
from §, their null inward normals are compatible with 
those used in the standard analysis. 

In this section, the subscript § is used to indicate the 
value of a quantity at § rather than at an arbitrary point 
of 3\f. However, having found the expressions for all quan- 
tities of interest in terms of data at S, in later sections 
almost all computations will be expressed in terms of 
these data, and the subscript will be omitted. 



parallel along l a mean the spin-coefficients e, k, 7' and r' 
all vanish. The choices so far apply to the spin-frame on 
INT, but it will also be convenient to restrict its behavior 
to first order off [NT. We require that V[ a Z b ] = 0, which 
implies r = a + ft. 

The restrictions e = 0, 7' = 0, and r = a + ft break 
the strict boost invariance of the Geroch-Held-Penrose 
calculus. However, a modified invariance still holds. If 
we consider rescaling o A H> Xo A . i A h-> A -1 /," 4 , where A is 
a non-zero complex- valued function on DsT which depends 
on the generator only, then the conditions e = 0, 7' = 
are preserved. If we consider that, accompanying any 
such rescaling the spinor field o A is changed to first order 
off N by D'o A h-> X~\D'o A - X~ 2 (BXX)i A ), then we find 

t ^ AA _1 (r+(AA)- 1 gAA), a+ft i-> AA _1 (a+/3)+A~ 2 3AA 
and t = a + ft is preserved. Since this modified trans- 
formation law for the spinor dyad differs from the simple 
rescaling of the GHP scheme only by altering the deriva- 
tive D'o A by a multiple of l a , the only spin-coefficient in- 
homogeneously affected on N is r; also the operators 3, 3' 
on N retain their usual GHP transformation rules. (The 
behavior of t adopted here is natural within the context 
of the characteristic initial- value problem; cf. (Tfl. 

Let s be an affine parameter along the outgoing null 
geodesies normalized so that l a \7 a s = I and vanishing 
at S. We may think of s having Geroch-Held-Penrose 
type {p,q} = {—1,-1} (given our restrictions Do A = 
0, Dl a — 0). (It is conventional to use r for an affine 
parameter if the boost freedom in the spin-frame is fixed 
to give the standard Bondi-Sachs asymptotics. However, 
as we shall have not fixed the spin frame in this way, we 
use s to avoid potential confusions.) 

All of the computations are built on the integration of 
the optical equations, which can be written as 



D 



p 


a 




P 


a 


a 


P _ 




a 


P _ 



(3) 



since v&n = and $nn = 0. The solution to this is 



A. The Integration Scheme 

Let S be a spacelike two-surface of spherical topology 
in a vacuum region of space-time, and let o A , i A and 
l AA ' = o A o A ' , m AA ' = o a l a ' , m AA = l a o a ' , n AA ' — 
l a l a be a spinor dyad and vector tetrad associated with 
it, so l a is the outgoing null congruence and n a is the 
ingoing null congruence. Assume next that this surface 
were embedded in a space-time with the same first and 
second fundamental forms, and with the same values of 
*f> n for I < n < 4 at §. (Actually, the value of ^4 will 
not enter.) We further assume that o A , l a are propagated 
parallel along the outgoing null congruence from S, and 
that ^>o vanishes along this congruence. This leaves a 
freedom o A n> Xo A , l a \-> X~ 1 l a where A is a function of 
the generator only. 

We have then that p = ~p on the null congruence, and 
also that p' s — p' s , where the subscript indicates restric- 
tion to §. The conditions for o A 7 l a to be propagated 



P ° 
a p 



1 
I 



P ° 
a p 



(4) 



Notice that these matrices, for different values of s, all 
commute, since they lie in the commutative algebra gen- 

" a~ 



erated by the identity and the single matrix 



a 



It 



is this which allows the integration of the system explic- 
itly. 

It is easy to see that p and a will be non-singular for 
s > iff p§ < — |<r§| < 0, which is the condition that no 
conjugate points develop. In practice one expects ps to 
fall off as 1/r and as to fall off as I /r 2 where peeling holds 
(where r is the affine parameter in the physical space- 
time, not the mathematical stand-in with ^0 = 0), so for 
large enough surfaces § which are close enough to spheres 
one should have |<7g| <C \ps\- This means that for such 
surfaces it should be a good approximation to neglect the 
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effects of a relative to p on propagation outwards. In this 
approximation, we have 



p 


a 




a 


P _ 





p§(l-Sp§) 1 <7§(l-Sp S ) 2 

ag(l - sp s y 2 pg(l - spg)- 1 



(5) 



Ps(l- 

to indicate the con- 



Here and throughout, we shall use 
dition | er s | "C |p§|. 

The remaining equations are integrated successively, as 
follows. (The results are given below, Tables HI HO ) From 
the equation 



-D*l = 4/9*1 



(G) 



one finds *i. With that in hand, one takes up that for 
the spin-coefficients a and /3, which we write as 
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a 




P o 




a 
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J_ 




a P _ 




.P. 


+ 





(7) 



This, and others to follow, can be integrated using the 
result 



exp 



f 


p ° 


(«') da' = 


P G 




p a 


10 


a p _ 


° P _ 




a P 



(8) 



(Recall that the matrices in the exponential all commute, 
so there is no need to take a path-ordered exponential.) 

One next takes up the transport of the operators d, 
5' up the generators. Lie transport along l a , the vector 
tangent to these, establishes a canonical diffeomorphism 
of the outgoing null surface INT with §x{s|s>0}. Using 
this diffeomorphism, we may extend 6§ = trig's/ a from its 
definition on § to N; this is equivalent to extending it 
by requiring it to be Lie transported along l a . If we put 
TOg = Am" + Bm a + Cl a , then using the standard spin- 
coefficient commutators we find the conditions for being 
Lie-transported are 





'A' 




P 




a " 




" A ' 
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B 




a 




p 




B 


, (9) 




C 




.-03 + 


a) - 


08 + a) 




C 




which can 


be integrated 


using 


J8]) and the initial condi- 


tions A = 


1, B = 


0, c = 


0. 









With the results of these integrations, one can find the 
remaining quantities of interest. For *2, we integrate 



L>* 2 = (6' - 2a)*! + 3p* 2 



(10) 



With this, we can integrate the transport equation for 
the optical scalars of the ingoing congruence: 
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(11) 



In fact, of these, we only need a' as a datum for the next 
equation, and we need to note that pf and a 1 will have, 
by virtue of ©, the asymptotic behavior 0(s~ 1 ). Then 
one can integrate 



L>* 3 = 5'* 2 + 2,0*3 + 2cr'*i . 



(12) 



Recall that we are interested in the asymptotic forms of 
those spin-coefficients and operators necessary to define 
the asymptotic spinors and twistors, and also the curva- 
ture quantities *i, *2, *3. Table fl] gives these results. 
The asymptotic values of the quantities a, /3, A, B, C 
enter into the determination of the asymptotic spinors 
and twistors, and there is little point in discussing them 
before those results are at hand. The asymptotic values 
of *i, *2, *3 are of some interest, however. For *2 and 
*3 , we see a leading term which is the value on § multi- 
plied by the power of (1 — sp) dictated by peeling, plus 
correction terms. One can check that, if the usual peel- 
ing assumptions hold in the physical space-time, then 
the scaling of these correction terms with the position of 
S is subdominant to that of the leading term. Thus the 
terms represent finite-size effects due to the distance of S 
from 3 + , and, despite their complexity, should be small 
for large enough §. The complicated nonlinear forms of 
these corrections are also of some interest. 

While we expect the results of Table U to suffice for 
most numerical work, it is of some conceptual interest 
to extend these by dropping the assumption |cr§| <C |p§|- 
This is done for those quantities needed to determine 
the asymptotic spinors and twistors in Table [TTJ Expres- 
sions for the curvature quantities *2 and *3 in this case 
are prohibitively lengthy; if any are needed, the best ap- 
proach is to do the corresponding integrals as power series 
in as/ ps, keeping as many terms as required. 



B. Transformation to a Newman Unti Frame 

The asymptotics of the tetrad and curvature compo- 
nents computed above can all be examined as s — > +oo, 
and compared with the requirements for a Bondi-Sachs 
frame, as given for example in ref. [8j. There are three 
sorts of adaptations which are necessary to make the 
tetrad accord with the standard formulas. 

One of these is to shift the zero of s to eliminate the 
s~ 2 term in the asymptotic expansion of p. However, this 
only alters subdominant terms in the expansions, and we 
will only need the dominant terms, so we omit this. The 
second change is to replace the spinor i A by one which 
becomes tangent to future null infinity as s — > oo; this 
can be accomplished by a null rotation. The final change 
is to break the local boost invariance of the dyad so as 
to achieve the standard asymptotic scaling p' ~ (2s) -1 . 
This requires solving an elliptic equation on § (equivalent 
to conformally uniformizing the sphere). 

Because of the additional computational resources re- 
quired to solve the elliptic equation on the sphere, we 
will distinguish here between the two stages of the pas- 
sage to the Bondi-Sachs frame. When the frame as been 
adjusted by a null rotation so one spinor is tangent to 
null infinity, we call it a Newman-Unti frame and in- 
dicate it by the postscript NU; after the rescaling, the 
Bondi-Sachs frame is indicated by B. Formulas valid in 
Newman-Unti frames will thus automatically be valid in 



s 



TABLE I: Asymptotic forms of the relevant spin-coefficients, coefficients of operators, and curvature, in the parallel-transported 
frame, in the case |os| -C \ps\- The left-hand column gives the asymptotic form of each quantity X for large affine parameter s 
in terms of a leading coefficient X°; the coefficient is given in the right-hand column. In the right-hand column, all quantities 
are to be evaluated at § (and so would, in the notation of this Section, ordinarily carry the subscript S, but this is omitted 
here) . 

Quantity X Leading coefficient X° , assuming \crg\ -C ps| 

A ~ A°s -p 

B~B°s 

C~C°s a + /3 

p~p°s _1 -1 

0-2 / \-2 

a ~ a s \~P) a 

(-p)- 4 *i 

* 2 ~ *as" 3 (-p)- 3 {* 2 -p _1 3'*i + 2p- 2 *i9'p + 2(3p)- 2 |*i| 2 } 

(-p) _2 {*3 + P _ V'*i + (l/6)g'[-4p- 3 *i9'p + 3p" 2 g'*i -p~ 3 |*i| 2 -p _1 *2] 
*3~*is~ 2 + (30p 4 )- 1 (g'p)[15*ig'p-10p9'*i+4|4'i| 2 + 15p 2 *2] 

+ (30p 4 ) _1 [15|*ii 2 g'p + 4\E' 1 ('*T) 2 - 10p*ig'*i + 15p 2 *T>f2]} 



TABLE II: Exact leading spin-coefficients and coefficients of operators in the parallel-transported frame. In the right-hand 
column, all quantities are to be evaluated at the surface S (and so would, in the notation of this Section, ordinarily carry the 
subscript S, but this is omitted here). 
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■ log 
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Bondi-Sachs frames. We shall give the transformation to 
a Newman-Unti frame in this section; the next section 
will cover Bondi-Sachs frames. 

The spinor l a representing the ingoing null congruence 
has been fixed by the geometry of the two-surface §, and 
propagated outwards along null geodesies; there is no 
reason to expect it to be asymptotically tangent to null 
infinity. We must therefore anticipate making a null ro- 
tation i A — »■ t^u = t A + Q° A in order to achieve this. To 
preserve the parallel-propagation condition Dl^ v = 0, 
we shall need DQ = 0. 

The equation determining Q comes from the fact that 
(with the correct choice of t^y) the quantity r must van- 
ish at least as 0(s~ 2 ) as s — 5- oo (cf. UGl). ( Note that 

this means Q will not transform homogeneously under 



rescalings of the dyad.) Since we have r = o a D'oa, we 
find t nu = o A (D + Q5 + QS')oa = r + Qa + Qp. Setting 
the lead asymptotic term of this to zero and using the 
asymptotic forms of p, a and r, we see that 



i- L J —a — p 



Now let us work out the operator S for this frame, 
which we denote <5nu = 6 + QD. We first invert the 
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definition of 5s to get 



III. BONDI SACHS FRAMES 



= (i^i 2 - i5i 2 r ] 

~ -s-\pl-\a s 
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-B 
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-CD 
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P -a 




"-5s 


- CD 


. ~ W p 


s 




-CD 



(14) 



where 


~ denotes asymptotic equality as s 


->• • 


foo), 


so 


£>NU 


~ -s-\pl -Wsl 2 )- 1 


P -a 
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CD 
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D. 



Inserting the asymptotic values from Table |TT1 we find 
that the coefficients of D cancel, and we are left simply 
with 



ft 



NU 



\pI-\°s\ 2 Y 
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P -o 








. ~° p 
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This cancellation is not a surprise; it arises from the fact 
that the s = const surfaces can be taken be the ones 
approaching cuts of 3 + . 

We now find the spin-coefficients qnu, /3nu m the 
Newman-Unti frame. (We shall not use these explicitly, 
but we give them both for completeness and for pur- 
poses of comparison with other work.) The transforma- 
tion rules for these are 



In this section, we complete the transformation to a 
Bondi-Sachs frame, and also establish some of the calcu- 
lus of these frames which will be required for the analy- 
sis of energy-momentum and angular momentum. The 
Bondi-Sachs frames are essentially equivalent to the no- 
tion of an "asymptotic laboratory frame." 

From this point on, we will omit the subscript § for 
spin-coefficient quantities (including operators 8, 5', 9, 
3') at the two-surface, unless explicitly indicated other- 
wise. Quantities considered at other points on 3Sf will ei- 
ther be in the Newman-Unti or the Bondi-Sachs frames 
and will be indicated by subscripts NU or B. 



(15) 

A. Transformation to a Bondi-Sachs Frame 

We recall that the Newman-Unti frame established in 
the previous sections is very nearly a Bondi-Sachs frame; 
what remains is to adjust the spin frame, or equivalently 
conformally transform, in order to make the s = const 
cross-sections unit spheres. 

Let us begin with the metric structure. This is char- 
acterized by the intrinsic 5-operators of the s = const 
surfaces. Equation (|16|) expressed these in terms of the 
structure at §; let us put <5nu = M a V a , so 



Q!NU 

/3nu 



a + Qp 



(17) 
(18) 



Using the results from Table |TT1 we find 



-^-ksfrH-piS+^+a/a)*!} s (19) 



v: 



§ -\<J^)- l {pP-aa-{l/2)^x} 



(20) 



results in the approximation |cr§| <C \p%\ are listed in 
Table HU 

The fact that Q remains bounded as s — > oo (in 
fact, Q is constant along the generators of 3Sf), together 
with the peeling property, means that transformation to 
the Newman-Unti frame does not affect the asymptotic 
forms of \&xi ^2, ^3- Also pnu = P an d ontj = fT ' smce 
these only involve derivatives of o A . 

The transformations for the optical coefficients for the 
ingoing congruence are more complex. We have p'^u = 
ff - 3Q - Q 2 <t, a' m = a' - d'Q - Q 2 p. We shall not 
need these, but for completeness their asymptotic forms 
are given in Table Mil 



M a 



-s^ip 2 - \a\ 2 )- 1 (pm a - afn a ) . 



(21) 



Then -2M< Q Af ; gives the intrinsic inverse metric on the 



-2M (a M b) , 



s = const surfaces. The intrinsic metric is 
where M a , M a are not defined via lowering with g a b, 
but via the duality relations Af a M° = -1, M a M a = 0; 
explicitly 



s(-pm a - am a ) 



(22) 



and 



2M (a M b) 



-2s 2 (pm (a + o-m {a )(pm b) + am h )) . (23) 



Under a change of scale o A \-> Xo A , l a m> X^ 1 ^, The 
factor (pm^ a + amu){pfn, b ) + cfmf,)) will be multiplied 
by |A| 4 ; we may therefore regard the change of scale as 
equivalent to a conformal transformation by 51 2 = |A| 4 . 
More precisely, for each choice of scale we have a family 
of metrics on the s = const surfaces which are scalar 
multiples of each other; when we change the scale we 
also change the choices of the s = const surfaces (since 
s 1 y I A I 2 -s) , and it is the metrics on these (pulled back 
to S) which are conformally rescaled by \X\ 4 (relative the 
pull-back of the metric on the original surface at the same 
numerical value of s). 

The conformal structure is characterized by the com- 
plex structure. We may introduce a complex stereo- 
graphic (antiholomorphic) coordinate [l8[ £ on the s = 
const surfaces by requiring Af a V a £ = 0, that is, 



(fiS - o-5')( = 



(24) 
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TABLE III: Leading forms of the spin-coefficients and operators in the Newman-Unti frame, under the assumption |<rs| -C |ps|- 
Each of these quantities falls off as s -1 with the corresponding coefficient. Thus qnu ~ s _1 «°, etc. In the right-hand column, 
all quantities are to be evaluated at S. 

Quantity Value 
S -p _1 5s 
S' -p~ 1 8' s 
oP -p- 1 {-p + (2p)- 1 *~i} 

P° -p- 1 {/3-(2p)' 1 *i} 

p" 1 { - p' - 2(3p 2 )- 1 |*i| 2 + (2p 2 )- 1 (9'* 1 +9*7) - (2p)" 1 (*2 + W 2 ) - 5(6p 3 )- 1 (*i9'p + *T3p) 
(p')° - (2p)- 1 (g- + g'r) + (2p 2 )- 1 (-3p + rg'p) + (2p 2 )- 1 (* 1 r + *Tr) + ^p 2 )" 1 (r9a + rOV) 

- (6p 3 ) _1 (*iga + *79'cr)} 

(<x')° -p- 1 {a' - 3((2p 2 )- 1 *7 - p-V) - P((2p 2 )- 1 *T - p-V) 2 } 



and requiring that £ be regular over § except for a simple 
pole (equivalently, that vanishes at a single point and 
in the limit of approach to this point its argument has 
winding number —1, the minus sign on account of its 
antiholomorphic character) . 

The coordinate £ is unique up to a fractional linear 
transformation. In order to keep its interpretation as di- 
rect as possible, when § is approximately a round sphere 
the coordinate £ should be taken to be close to a stere- 
ographic coordinate cot(6*/2) on S. One way to fix 
the freedom would be to require the pole to lie on the 
+z coordinate axis, the zero to lie on the — z axis, and 
the point £ = 1 to lie on the +x axis. With these choices 
we effectively fix an asymptotic "laboratory frame." (The 
time axis is fixed by the requirement that |£| = 1 be a 
great circle.) 

For numerical work, it may be more convenient to re- 
cast Eq. (|24l) in terms of regular quantities. If we let 
Co be any smooth function on § with a simple pole and 
simple zero of the required type and write £ = ££o, then 
£ is smooth over § and satisfies the everywhere-regular 
equation M a V a ( = -(M Q V a log Co- There will be 
a one-complex-dimensional space of everywhere-regular 
nowhere-vanishing solutions to the equation for £. These 
solutions will give £ = £Co the same pole and zero as 
Co; thus, if these have been chosen as in the previous 
paragraph, one has simply to adjust the multiplicative 
constant in £ to achieve the final normalization £ = 1 
on the +x axis. We shall assume from now on that a 
solution C to (fM| has been found. 

Now let us turn to the metric structure. As is con- 
ventional, put M a = — P d(, so that — 1 = M M a ~ 

P'V 1 ^ 2 - H 2 )- 1 ^' ~aS)C ~ P~ 1 a- 1 p-H'Z and 

P ~ -s-V^C- (25) 

Then the metric is ~ -2\P\^ 2 dC,dC . 

Now let us consider a change of scale to achieve a 
Bondi-Sachs frame. Let this be o A o4 = Xo A . We 



may keep £ unchanged (it is a conformal invariant); we 

have then P m> Pb = — «b Vb^MbC ~ —A s Jj 1 p~ 1 5(. 
We rescale the metric to a sphere of radius sb with 
P B = 2 -1 / 2 s B 1 (l + |C| 2 ). This will align the time axis 
of the Bondi-Sachs system with that of the laboratory 
frame, and give the standard spin frame adapted to the 
coordinate £. We then find 

A 2 = -2 1 / 2 p- 1 (l + |C| 2 )" 1 <5'C. (26) 

With A known, we may read off the value of any spin- 
and boost-weighted quantity in the Bondi-Sachs frame 
from its value in a Newman-Unti frame. (Note that the 
original o , l a may be any spin frame for which m a is 
tangent to §; Eq. (|26j) provides the correct transformation 
to the frame adapted to £.) 

In what follows, we shall need the shear in a Bondi- 
Sachs frame. We have, from Table UIl that ctnu ~ 
s~ 2 (p 2 — \a\ 2 )~ 1 a. Inserting the appropriate rescalings, 
we find 

vb ~ s B 2 X%- 3 (p 2 - |a| 2 )-V (27) 

with A given by (|26|l. In particular, the Bondi shear is 
the coefficient of Sg 2 , that is 

£r B = AA 3 (p 2 — | cr | 2 ) 1 cr . (28) 

Finally, we remark that an alternative (and somewhat 
more traditional) route to fixing the conformal factor is to 
require that the Gaussian curvature of the s = const sur- 
faces becomes asymptotically constant. Since the Gaus- 
sian curvature is ~ — 2/9nuPnu anc ^ '° NU ~ — s_1 ; this 
leads to the requirement that p^u ~ (2s) -1 . One can use 
Table Mil and the transformation rules discussed at the 
beginning of section Hi Al to write this as a second-order 
partial differential equation for |A|, which is equivalent to 
the usual equation for finding a conformal transformation 
uniformizing the Gaussian curvature. The formulas for 
this are rather more complicated than those given in the 
present subsection, however. 
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B. The Angular Potential for the Shear 

The Bondi shear, being a spin-weight two quantity, 
admits an angular potential Ab such that 

43|A B = <r° (29) 

(or equivalently 9 b Ab = &b)- The use of this potential 
facilitates the computation of the Bondi-Sachs energy, 
and the potential also plays a central role in the analysis 
of angular momentum. The electric and magnetic parts 
of the shear are a c \ — 8|5RAb and <7 mag = z3 b 3Ab- 

Equation ((29]) is easily solved when Pb = 2~ 1/ ' 2 s B 1 (l + 
|£| 2 ). In this case, it can be written as 

(1/2)^(1 + CC) 2 ^A B - a° , (30) 

and a Green's function for the operator can easily be 
derived from the relation d^C — (£, C) (where the 
right-hand side is the usual 6- function in the £-plane). 
We find 

Ab(C,C) = J G((X;(,t)<JB(Lt)d<;Adt/(2i), (31) 
where 

G (C,c i c,r) = -.-K-0-'( T ^ ?F - I ^) 

(32) 

and 

(2i)- 1 d(Ad(= (1 - \<j/p\ 2 )\6(\ 2 dS. (33) 

An alternative means of solving Eq. (|29|) would 
be to resolve <r B into spin-weighted spherical har- 
monics 2Yj,m and then use the relation 9| oYj. m — 
(l/2)y/(j - + + 2) 2 Yj, m to infer the corre- 
sponding resolution of Ab- Thus one would have 

4m = 2 1/2 j2Y~ (l + ld 2 )- 1 ^^^^, (34) 

X j>m = 2((j - l)j(j + l)(j + 2))- 1 /V° m (j > 2) , (35) 
Ab = Xj, m oYj,m ■ (36) 

The terms A J)TO with j = 0, 1 are freely specifiable. We 
take these terms to vanish, which will simplify the coor- 
dinatization of the twistor space, below. 

The phases of the spin-weighted spherical harmon- 
ics depend on the spin-frame. There are two common 
choices: that adapted to the complex coordinate and 
that adapted to 0, <f>. Because the analysis has been given 
here in terms of it is that spin-frame and those spher- 
ical harmonics which are used in Eq. (|34]) . To use the 
harmonics with respect to 9, <f> one must, besides replac- 
ing _ 2 !j, m (C,C) with -2Yj : m(0,(l)), also include a factor 
of (C/C) 2 in the integrand of (fBTI) . (Cf. the appendix; see 
ref. [14[ for a detailed discussion of the harmonics.) 



IV. ASYMPTOTIC TWISTORS AND SPINORS 

The Bondi-Sachs energy-momentum is a covector in 
a certain vector space, the space of asymptotically con- 
stant covectors. This space is most easily constructed 
from the space of asymptotically constant spinors. Simi- 
larly, the twistorial angular momentum is defined on the 
space of asymptotic twistors. In fact, the asymptotic 
spinors are naturally defined in terms of a canonical fi- 
bration of twistor space, so we shall start by construct- 
ing the twistors and then specialize to the spinors. We 
conclude this section by giving the kinematic twistor, in 
terms of which the energy-momentum and angular mo- 
mentum will be defined. 



A. The Twistor Space 

The twistor space T(C(S)) of the cut C(S) of null in- 
finity associated with § is the set of solutions of the two- 
surfaces twistor equation at C(§). These equations are 

3b^b - (37) 
sSb^b = (38) 

where cj b = s _1 oj° is rescaled to attain a finite limit at 
C(S), and s9 B tends to an operator depending on angle 
only; cf. Eq. ([25]) . There is a four-complex-dimensional 
space of solutions to these which we shall give shortly. 
For completeness, we note the forms of these equations 
in terms of the spin-coefficients at § are 

- {p& -m+ (l/2)¥0w° = (39) 
~(pg-o-9' + (l/2)*i)w 1 = (7w°. (40) 

(The minus signs are included because p is negative.) 

B. Solving the Twistor Equation 

Solutions to the twistor equation are easily found. The 
equation (|37|) for ui° has as its space of solutions the 
spherical harmonics of spin- weight —1/2; thus 

^ B = 2 i/2 (i + ici 2 r 1/2 (z 3 + ^ 2 c), (4i) 

where Z 2 and Z 3 are constants. To solve the remaining 
equation, we adopt a device of K. P. Tod and set 

w B = w B 3 B A B - A B 9 B w B + £b , (42) 

where Ab is an angular potential for the shear (subsec- 
tion |niBj, and £b is a spin- weight +1/2 quantity to be 
determined. We note that 

Qb^b = ^ + \C\ 2 r 1/2 (Z 2 -Z 3 (). (43) 

Then the remaining equation (|38l) is equivalent to 9 b £b = 
0, and the solutions to this are 

iB = -i{l + \C\ 2 r 1/2 {-Z Q + Z 1 Q. (44) 
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Thus (Z°, Z 1 , Z 2 , Z 3 ) coordinatize the twistor space; the 
factors have been chosen to make them accord with those 
induced from the standard Cartesian basis of Minkowski 
space if the cut of null infinity is got from the light-cone 
of the origin (cf. [lij], section 4.15 and Q, section 6.1)[l9|. 



C. Structures on Twistor Space 

There are three important structures on twistor space: 
a fibration, an infinity twistor, and a certain reality struc- 
ture. The fibration and the infinity twistor allow the 
definition of asymptotic spinors; the reality structure de- 
fines the null geodesies which play the roles of origins for 
the definition of angular momentum. Finally, the real- 
ity structure and the infinity twistor combine to define a 
certain twistor operation, the "hook," which enters in the 
definition of energy-momentum. 



1. The Fibration, the Infinity Twistor and Spinors 

The most primitive and important structure on twistor 
space is the fibration, which is defined by simply keeping 
the ui° field of the twistor. This is just (Z°, Z 1 , Z 2 , Z 3 ) h-> 
(Z 2 ,Z 3 ) in our coordinates, since w° is specified by Z 2 
and Z 3 . We see then that the space of fibers is a two- 
complex-dimensional space; it is naturally identifiable 
with the space of (dual, primed) asymptotically constant 
spinors E>a>- Asymptotic spinors of other valences and 
asymptotic vectors are tensors are defined as usual by 
tensor operations from §a' ■ 

Closely related to the fibration is the infinity twistor 
I(Z,Z) = I al3 Z a ZP. In our coordinates, it is simply 
given by 



I(Z,Z) = Z 2 Z 3 - Z 3 Z 



73 62 



(45) 



This evidently defines a skew form on Sa>, which rep- 
resents the asymptotically constant spinor e 



A'l 



Its 



negative inverse is ca'B' 4 , the spinor €ab^A'B' represents 
the asymptotic metric g a b — 9AA'BB'- One often puts 
Z 2 = irni . Z 3 — 7Ti / and then one has 



with e° 1 = 1 as usual. 

One can also define spin space E> A directly in twistor 
terms, as the kernel of the fibration; thus E> A is iden- 
tified with the set of spinors whose coordinates are 
(Z° , Z 1 ,0,0). However, it will be more natural for us 
to work with Sa 1 , especially when we take up evolution. 



I{Z,Z) 



C A'B' - 

e tta'TTb' 



(46) 



2. The Reality Structure 

The twistor space for Minkowski space is equipped 
with a sesquilinear form of signature + + — — whose 



zero set is the set of real twistors. In general relativity, 
there is also a reality structure, but it is more nonlinear 
when spin is present. The analytic manifestation of this 
is that the candidate expression for the norm 



— o' 

- U! 



B"B^B 



J B U B^B 



Q" B ) (47) 



is not in general constant over S. It turns out that we get 
a good theory of angular momentum, however, by simply 
evaluating this expression at the point on the sphere for 
which Ug, or equivalently vanishes. (There always 
will be a unique such point, unless the field w° vanishes 
identically, in which case the result is taken to be zero 
by a continuity argument.) We will denote the point at 
which u>° vanishes by 7(w°); the restriction of (|4"Tf to 
j(u>°) by $(Z). A twistor is real iff $(Z) = 0. The 
real twistors correspond to real null geodesies meeting 
null infinity; they take the place of space-time points as 
"origins" for the evaluation of angular momentum. 
Using the formulas above we find that 



$(Z) = -23A|S B £ B | 2 + K&b^ + £o>b': 



7 (iD°) 



(48) 

The stereographic coordinate of 7(a) ) is ( = —Z 3 /Z 2 
and a brief calculation gives 



Using these, we find 



$(Z) 



-23A 



\Z 



2 1 2 



\z 



3|2n 



-Z Q Z 2 + Z 1 Z 3 + Z 2 Z° + Z 3 Z X 



(50) 



Here the second line would be the usual twistor norm in 
special relativity; the contribution on the first line is an 
essentially general-relativistic effect. We recall that 3A is 
the angular potential for the magnetic part of the shear; 
it is this magnetic shear which distorts the twistor norm. 
(The magnetic shear is found to be the j > 2 components 
of the specific — that is, per unit mass — spin.) 



3. The Hook Operation 

While, as discussed above, the special-relativistic 
twistor norm does not extend to general relativity, 
enough of the structure does survive that a certain an- 
tilinear operation does carry over to general relativity. 
This is the hook of a twistor Z a , denoted by Zp . (Ow- 
ing to the non-existence of a norm, the quantity Zp does 
not have any separate meaning for us, but we shall keep 
the special-relativistic notation.) If we set Z a — I^ a Zp, 
then the definition of the hook is 



-0 

= 

- 1 

U) = — 



—0' 

UjJ 



(51) 
(52) 
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In our coordinates, the hook operation is 
(Z°,Z\Z 2 ,Z 3 ) H> (-Z 3 ,Z 2 , 0,0). (The hook of 
any twistor is an element of E> A .) 



D. The Kinematic Twistor 

In twistor theory, the energy-momentum and angular 
momentum are encoded in a kinematic twistor A(Z) = 
A aP Z a Z^ defined by 

A(Z) = - i (4 7 rG)- 1 j{^{u%) 2 + 2^ u ^nu 

+ *f J (^ u ) 2 } d5NUj 

(53) 

where the limit as the surface tends to null infinity is un- 
derstood. The kinematic twistor satisfies the Hermiticity 

property A a/3 Z T^ 1 = A a/3 Z a I'<^Z 7 , which in our co- 
ordinates is 

A 00 = Am = A u = (54) 
A 02 = —Ai3, A u = A^, A 03 =Aq3. (55) 

To work out the kinematic twistor explicitly, let us first 
insert the asymptotic forms of the quantities; we find 



(56) 



of asymptotically constant spinors. In twistor terms, this 
is 



A{Z) = -i{AKG)j{p 2 -\o-\ 2 ) 

x {*?(£ ) 2 + 2*°cA/ + ^(w 1 ) 2 } dS , 



in terms of data on S. 

We next express the twistors in terms of their forms in 
the Bondi-Sachs frame, taking into account the rescaling 
relative to the Newman-Unti one (recall that oj°, to 1 have 
Newman-Penrose types {0,1}, {1,0} respectively): 



A(Z) = -»(47rG) I 2 {Col) 



+2^° 2 {XXr 1 Cj° B Lji + *°A- 2 (c4) 2 } (p 2 - \a\ 2 ) dS . 

(57) 

At this point, we may substitute the explicit forms of 
the solutions of the twistor equation given in this section 
to compute the components A a p of the kinematic twistor. 
Because of round-off errors (and also the approximation 
rr | <C \p\, if used), the numerical computation of A a p 
directly from (|57|) could fail to satisfy the Hermiticity 
conditions (j54")h (155)) ; we therefore enforce these condi- 
tions at the levels of the integrands. The results of this 
are given in Table IIVI 



ENERGY-MOMENTUM 



The Bondi-Sachs energy-momentum P a is most natu- 
rally viewed as a function P AA ttatta' on the space E>a' 



P TT A TT A , 



A al3 z a r?z 1 , 



(58) 



where the right-hand side represents the contraction of 
the kinematic twistor once with Z a and once with its 
hook 7 7 ^Z 7 ; the result (|58| is real depends only on the 
projection 7r A' of the twistor Z; that is, if the twistor 
components are (Z°,Z 1 ,.Z' 2 = ttq',Z 3 = ny), the choice 
of Z° and Z 1 is immaterial. 

Using the formula for the hook map in coordinates 
(just below Eq. (f5"2"j) ). we find explicitly 



(59) 



pOO' 


pOl' 


= 2- 1 ' 2 


pt 4 


_pz 


px 


+ ipy 


plQ' 


P w 




px _ 


ipy 


P t 


- p z 



Al2 
-Aq 2 



Ax3 
-Aw 



VI. ANGULAR MOMENTUM 

In special relativity, the angular momentum is a 
position-dependent skew tensor field M a b(x), or equiva- 
lently a spinor p A 'B'(x), where M ab (x) = PA'B'(x)e A B + 
~Pab(x)£A'B' ■ In general relativity, as is well known, 
there is no good set of "asymptotic origins" for the mea- 
surement of angular momentum, and thus the special- 
relativistic treatment does not apply. 

A good treatment of angular momentum in general 
relativity is possible, if we adjust our perspective a bit. 
We first note that if 7 is a null geodesic in Minkowski 
space and tta 1 is a tangent spinor to it, then the compo- 
nent p, A B (x)tta'T^b' of the angular momentum is (by the 
change-of-origin formula) independent of long as x 
lies on 7. Thus we may think of the angular momentum 
either as a spinor-valued function on space-time, or as 
p(l,^A') = p A B (x^a'TTb' j a function on the space of 
null geodesies together with their tangent spinors. While 
angular momentum does not extend to general relativity 
as a spinor-valued function, it does extend as a function 
of the null geodesic and the tangent spinor. Indeed, the 
expression is very simple: we have 



ii(rr,T A >) = (2i)- 1 A aP z a z 13 . 



(60) 



where Z <->■ (7, tta') is the real twistor defined by the null 
geodesic 7 and the tangent spinor tta 1 ■ 

While the general-relativistic angular momentum is 
thus an extension of the special-relativistic concept, some 
features which are prominent on the special-relativistic 
case do not extend to general relativity, and others which 
are usually viewed as secondary become central in the 
general-relativistic setting. 

The root of this is that the general-relativistic an- 
gular momentum is defined on the space {(7,^/)} of 
null geodesies together with their tangent spinors. In 
the asymptotic regime, this space naturally fibers over 
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TABLE IV: Integrands for the computation of the potentially non-zero components A a p = Ap a of the kinematic twistor. Each 
term on the right is to be multiplied by the common factor — i(47r(?) -1 (l + |C| 2 ) _1 (p 2 — l°"| 2 ) an d integrated with respect to dS 
over the sphere. In this table, the Hermiticity conditions A\ 2 = Ai 2 , Aq 3 = A03, A02 ~ — Ai 3 have been enforced at the level 
of the integrands. 



Component Quantity to be multiplied by — i(4mG) 1 (1 + \C\ 2 ) 1 {p 2 — | cr | 2 ) and integrated with respect to dS 

A 02 = -M; 2 1 / 2 ^|A|- 2 K*§ + iQMpDcW^Xa) + 2- 1 i(? 2 - l)K(*gA- 2 A B ) - 2" 1 (f + 1)$}(*§A- 2 A B ) 

A 03 2^ 2 i\X\~ 2 ^° 2 + j5R(*iA- 2 (2 1 / 2 3BA B + <A B )) 

A 12 -2 1 / 2 l |C| 2 |A|- 2 ^ - i»(*gC(2 1/2 C3 B AB - A B )) 

A 22 _ _ 2*?r 2 f + 2 3 / 2 *§|Ar 2 C(2 1 / 2 C9BAB-AB) + ^A- 2 (2 1 / 2 C9BAB-A B f 

A 23 2*»A" 2 C + 2 1 / 2 ^|A|- 2 (2 3 '' 2 C3bAb + (Kj 2 - 1)A B ) + *iA- 2 (2 1 / 2 g B A B + CAb)(2 1 / 2 C9bA b - As) 

A 33 2*?A~ 2 + 2 3 / 2 *§|A|- 2 (2 1 /2g B A B + CA B ) + *gA- 2 (2 1 / 2 g B A B + CA B ) 2 



the space §^4/ of spinors, because there is a well-defined 
asymptotic spin space. This contrasts with the usual 
view of /i A B (x)tta'T^b' being defined on the spin bun- 
dle {(x, ft A')} of Minkowski space, where the base space 
is Minkowski space and the fibers are copies of Sa> . In 
practical terms, this means that while the component of 
the angular momentum /z(7, tta 1 ) m a direction corre- 
sponding to 7 and 7r a> will be well-defined, there will be 
no natural way of simultaneously varying 7 and tta' so 
the angular momentum is specified by a pure j — 1 rep- 
resentation [i A B of the Lorentz group; the angular mo- 
mentum will inevitably (if a magnetic part of the shear 
is present) have j > 2 parts as well. 

The reason for the appearance of these j > 2 com- 
ponents is that general relativity unifies the "ordinary" 
(J = 1) angular momentum with gravitational radiation. 
The j > 2 parts of the angular momentum correspond 
exactly to the shear (times the Bondi mass). Because 
there is no split of the angular momentum into j = 1 
and j > 2 parts with the appropriate geometric invari- 
ance (invariance under the Bondi-Metzner-Sachs group), 
one must, to get an invariant theory, consider all of the 
j > 1 parts of the angular momentum. 

While the angular momentum does depend on the 
pairs (7, tta')j the dependence is not arbitrary: there is 
a general-relativistic analog of the change-of-origin for- 
mula, which says that the angular momenta at different 
points in a fiber differ by appropriate multiples of the 
components of the energy-momentum. Thus the full in- 
formation in the angular momentum can be recovered 
by choosing any cross-section *y(tta') of the fibration and 
computing n(pf{KA')i'KA') as TA' varies. Because this is 
a homogeneous function (of degree two) in tta' , the es- 
sential information in the angular momentum is that in 
one spin- weight minus one function on the sphere. 

It is natural for us to choose the cross-section to be 
given by the congruence of null geodesies meeting the 
cut of null infinity orthogonally (that is, the congruence 
specified by l a ). This congruence then serves as a sort of 
origin for the computation. (However, the present pre- 
scription differs essentially from previous attempts to use 
cuts as origins.) Then the electric and magnetic parts of 
/i(7(7T J 4'), tta>) represent the energy moments and spa- 



tial angular momentum, respectively, of the system with 
respect to the asymptotic laboratory frame. 

We will also want a general-relativistic extension of 
the "polarized" form fi A B tta'^B'- This corresponds to a 
two-point function 

n((ry,ir A >),tf,*A>)) = (?{)- 1 A a f i Z a 2< i (61) 

on twistor space, where Z <-> (7,7^/), Z «-» (7, ~ka')- 
In special relativity we would have (2i)~ 1 A(Z,Z) = 
fi A B {x &v )iza i 'Kb' 1 where x av is any point on the world- 
line defined by "averaging" the geodesies 7, 7 with re- 
spect to the energy-momentum (x av = (x + x)/2 where 
x e 7, x e 7 satisfy (x a - x a )P a = 0) [9]. While in 
general relativity there is no similar invariant notion of 
averaging null geodesies, that is a limitation only on in- 
terpreting the origin of the angular momentum in direct 
space-time terms and not on its well-definition as a con- 
served quantity. 

If we fix a cross-section of twistor space, then the an- 
gular momentum /u((7(tta')) 7r A')i (tO^-A')) ^A')) can be 
thought of as a two-point function on the sphere. How- 
ever, the essential information in it corresponds to func- 
tions of one point on the sphere, not two. This is because 
the condition (f54"| implies that the higher-j terms enter 
only in tensor products with j = 1/2, s = —1/2 rep- 
resentations, that is, the angular momentum is in fact 
the symmetrized tensor product of a single spin-weight 
minus one-half function with an ordinary spinor. 

The intrinsic spin may be computed by passing to a 
boosted frame in which the time-axis is aligned with the 
Bondi-Sachs energy-momentum; then the magnetic part 
of the angular momentum is the spin. (The electric part 
of the angular momentum in this frame has a natural 
interpretation, too. The twistorial construction makes 
the cut appear as if it were a supertranslated cut in a 
stationary space-time; the electric part is the Bondi mass 
times this supertranslation.) 

A. Reporting the Angular Momentum 

Besides the need to accommodate j > 1 representa- 
tions, there is another issue to address in choosing how 
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to report the angular momentum, which is a trade-off 
between invariance and intuitive familiarity. This is- 
sue already occurs in special relativity, where the an- 
gular momentum is invariantly an element of the j = 1, 
s = — 1 representation, but when a reference frame is cho- 
sen we usually think of it as two spatial vectors (really 
the elements of a complex j = 1, s = representation). 
We shall opt for the familiar presentation, so that we 
can speak of the spatial angular momentum and energy- 
moment (both with contributions for j > 1) parts of the 
full angular momentum. 

Consider for a moment special relativity. Let t a be 
a unit future-pointing timelike vector, and let z a be a 
unit spacelike vector orthogonal to it. (For this para- 
graph only it will be convenient to regard z a as a vari- 
able direction on the sphere.) Then 2~ x / 2 {t a + z a ) is 
a future-directed null vector, say W A n A (normalized by 
t AA ,7r A TT A ' = 2-V2), and t/ A \ B '^ A = t A( - A 'z B ") B e AB . 
It follows that 

2 1 /^ A , Bl t A A 'jr A n B ' = M ah t^ A 'z B '^ B 

= M ab (l/2)(t a z»-(i/2)e abcd t c z d ) 
= (l/2)(M tz -iM xy ), (62) 

where (x, y, z) form a right-handed spatial triad. Thus, 
having fixed t a , as tt A i varies, one gets the energy- 
moment M tz and spatial angular momentum M xy in the 
direction z AA ' = 2 1 / 2 W A ir A ' - t AA ' it determines. 

We shall do the same thing in general relativity. We 
take for t a the time direction determined by the asymp- 
totic laboratory frame, we allow Z 2 = ttqi , Z 3 = iry to 
vary (normalized to | Z 2 \ 2 + \ Z 3 \ 2 = 1), and we take tta> — 
2 1 / 2 t AA 'if A , that is 7To' = 7fi, 7Ti' = —t^o- With these 
restrictions (i(ita') = ^{{l{^A>),^A'),{l{^A')^A')) be- 
comes a spin-weight zero function on the sphere, with 
fx + ~p, giving the energy-moment, and in — iJl giving the 
spatial angular momentum, in the direction determined 
by tta' (and t a ). These may be reported as real functions 
on the sphere, or resolved into spherical harmonics. 

(Of course, there is some freedom in choosing how to 
extend the terminology appropriate to a purely j = 1 
quantity to a j > 1 one. For instance, what one chooses 
to call the j > 2 energy-moments and spatial angular 
momenta could be taken to be some function of j times 
the ones used here. Such differences are unimportant 
here, since what we are interested in is simply extracting 
the invariant information.) 

B. Derivation of the Formula 

The twistorial formula for the angular momentum is 
simply 

H = (2i)- 1 A al3 Z a zP , (63) 



where A a p is the kinematic twistor introduced earlier and 
Z a , Z a are twistors whose null geodesies meet the cut of 
null infinity orthogonally (and satisfy \Z 2 \ 2 + \Z 3 \ 2 = 1, 
Z 2 =lP, Z 3 = -Y 2 ). The condition that Z a (say) meet 
the cut at a point is that the fields (cjg,Wg) vanish there; 
that the meeting be orthogonal means that the tangent 
spinor to the geodesic must lie in the o A / direction. 

Let the point in question on the cut have stereographic 
coordinate (. Then from the formulas (gTJ, ([32]). (03), 
(j44| for the twistor fields, we deduce that the conditions 
for 7 to meet the cut are 



C = -Z 3 /Z 2 , (64) 
X B (\Z 2 \ 2 + \Z 3 \ 2 ) = i(Z°~Z2 + Z 1 ^) . (65) 



The component of the tangent spinor in the l a direction, 
which we require to vanish, was computed in Eq. (29) of 
ref. [||. It is 



( | Z 2|2 + | Z 3| 2) -1/2| Z 2|-1 X 

{i(& B \ B )(\Z 2 \ 2 + \Z 3 \ 2 )Z 2 + 2- 1 / 2 {Z°Z 3 - Z 1 Z 2 )Z 2 ) , 

(66) 



where cJgAe is evaluated at (|64D . 

After a little algebra, we find the equations for Z°, Z 1 
in terms of Z 2 = ttqi , Z 3 = ny : 

Z° = -i(2 1 / 2 Z 2 Z3g^A B + |^ 2 | 2 A B )/Z2, (67) 
Z 1 = i(2 1/2 Z 2 d' B X B - Z 3 X B ) , (68) 

where A B and 9 B A B = 2- 1 / 2 (l + \(\ 2 )d\ B /d( are evalu- 
ated at £. Eqs. (|6"T)) . (f6"8")) determine the cross-section of 
twistor space as a bundle over spin space. 

The same analysis applies, of course, with Z a replaced 
by Z a . We note that ( = — 1/£ is the point antipodal to 
£ on the sphere, and the quantities A B , <3 B A B appearing 
in the formulas for Z°, Z 1 must be evaluated at this 
antipodal point. 

Substituting these formulas into Eq. (|63[) and collecting 
like terms, we find (with the normalization \Z 2 \ 2 +\Z 3 \ 2 = 
1) that 
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M = 2~ 1 (1 + Id 2 )" 1 { - [A 13 + (A u + A 03 )( + ^02C 2 ]2 1/2 3^A B (C, C) + [^03 + A 02 ( - A 13 ( - A 12 \{\ 2 ]\ B (C, C) 
+ [A02 - (A 03 + A 12 )( + A 13 C 2 ](C/C)2 1/2 g^A B (-l/C,-l/C) 

+ [A 12 + A02C ~ A 13 ( - ^ 03 |C| 2 ]A B (-1/C, -1/C) 
+iA 22 ( - iA 23 (\(\ 2 - 1) - iA 33 {} . 



(Here we have followed the standard convention of writ- 
ing a non- holomorphic function / of the complex coor- 
dinate C as /(£, £)•) While the formula is lengthy, this is 
mostly due to the appropriate factors of £, £ for weighting 
the components A a p of the kinematic twistor; the only 
complicated expressions are those involving Ab, which 
encode the shear. 

Thus Eq. (|69p allows one to report the angular momen- 
tum as a function of the direction (with 29fyt giving the 
energy-moment in the direction, and — 23/i the spatial 
angular momentum about the axis, specified by £). 

If the resolution of Ab in spherical harmonics is known, 
one can use it to find the resolution of /i into spherical 
harmonics, by identifying the explicit functions of £ in 
(|69f with particular spin- weighted harmonics and apply- 
ing tensor product formulas ("Clebsch-Gordan decom- 
positions"). The computation is lengthy but straightfor- 
ward using the formulas derived in the appendix; we find 

+ (2i) -1 y/2T/3(A 22 0Y1.1 + 2^ 2 A 23 Y lfi + A 33 Fi,_i) , 

(70) 

where the last three-terms are the AB-independent ones, 
and the coefficient jj,j tTn is a sum of terms, each of which 
is Xj, m : for f = j - + 1, m' = m - l,m,m + 1 
times a factor; these are given in Table [V] (for j even) 
and Table ED (for j odd). 

VII. EVOLUTION 

The usefulness of a definition of energy-momentum or 
angular momentum in radiation problems depends con- 
siderably on whether it admits a well-defined notion of 
evolution. At null infinity it is well-known that the 
Bondi-Sachs energy-momenta at two cuts can be com- 
pared. Many proposed definitions of angular momen- 
tum at null infinity took values in cut-dependent spaces, 
making tracking their evolution problematic; the twistor- 
based definition solves this problem. Here, however, be- 
cause we are dealing with energy-momentum and angu- 
lar momentum at large but finite spheres §, we must take 
up the problem anew, both for energy-momentum and 
angular momentum. 

Suppose we have a one-parameter family of surfaces 
§(rj) (for r\ in some interval J) foliating a timelike sur- 
face T, with r\ increasing towards the future (precisely, 



I 

we require v a \7 a r] > for every future-causal vector v a 
tangent to T). We may compute the energy-momentum 
and angular momentum on each of these; the difficulty is 
that these quantities are naturally defined on the null in- 
finities J + (rj) of the different auxiliary space-times M(r]) 
(each defined by taking = along the null hypersur- 
face N(?7) outwards from §(77)). In order to compare the 
energy-momenta and angular momenta for different rj, 
then, we must find a natural way of identifying the null 
infinities J + (fy) for different rj. 

We could express this as a problem of finding transition 
functions. The constructions above determine a preferred 
Bondi coordinate system (£,,, ( v , u v ) on U + (r?). (We have 
shown how to fix, for each §(77), a complex stereographic 
coordinate ( — ( n on J + (7/) , and we have chosen an asso- 
ciated Bondi-Sachs frame by fixing the factor Pb- This 
determines the Bondi retarded time u = up to a su- 
pertranslation; we fix this by requiring the preferred cut 
Q{rj) of J+(7 ? ) — the limit of }%) — have u v = 0.) Then 
we wish to find formulas for (C, )2 , £, )2 , u V2 ) in terms of 
(Ct?i : Cni ' u m ) ■ 

Structures Preserved by the Identifications 

I will discuss the way these identifications are made 
shortly. More important, though, is the question, What 
structures are preserved by the identifications? Were 
the surfaces §(77) actually cuts of the null infinity of the 
physical space-time, we should expect the usual struc- 
tures to be preserved and the identifications to be Bondi- 
Metzner-Sachs transformations. However, here we must 
expect a weaker structure due to finite-size effects. The 
extent to which these effects are significant should be 
interpreted as the extent to which the extraction sur- 
faces §(rj) are insufficiently distant to capture the full 
radiative structure (or, if the effects persists for arbitrar- 
ily large surfaces, the extent to which the Bondi-Sachs 
asymptotics fail for the physical space-time). Even with 
this interpretation, though, in order to precisely quantify 
the finite-size effects we must work out the structure in 
general. 

The usual intrinsic structure of null infinity may be 
regarded as determined by three elements: its set of gen- 
erators, each with an affine structure; a conformal struc- 
ture on the set of generators (it is a remarkable feature 
of the construction that the set does have a well-defined 
conformal structure); and the "strong conformal geome- 
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TABLE V: Terms contributing to the Ae-dependent part of the angular momentum proportional to oYj m for even j. Each 
term is the product of the component Xji m i in the left column by the quantity in the right column. (We understand the term 
is zero unless j' = 0, 1, 2, . . ., m' = — j', . . . ,j'.) 

Factor it multiplies 

2- 1 A )2 ((j-m + l)(j + m)) 1 /2 

-2- 1 (A 03 + A 12 )m 
-2-^13(0/ + m+l)(j-m)) 1 / a 



Component Xj 

Xj.m — l 

Xj J rn 
Xj,m+1 



TABLE VI: Terms contributing to the AB-dependent part of the angular momentum proportional to oY, m for odd j. Each term 

is the product of the component Ay >m / in the left column by the quantity in the right column. (We understand the term is 
zero unless j' = 0, 1, 2, . . ., m' = — j' , . . . , j' ■) 

Component Aj/ >m ' Factor it multiplies 

Aj-i, m -i 2-M 02 (i - 2)((j + m)(j -m+ l)/(4j 2 - l)) 1 / 2 

Aj+i,m-i 2-M 02 (j + l)((j - m + 2)(j - m + l)/(4(j + l) 2 - I)) 1 ' 2 

Xj,m 2 1 (Ao3 — A12) 

Ai-i, m+ i 2- 1 A ls (? - 2)((i - m)(j - m - l)/(4j 2 - l)) 1 / 2 

Ai+i,m+i 2-M 13 (j + l)((j + m + 2)(j +m+ l)/(4(j + l) 2 - l)) 1 / 2 



try" (which links the scales of vectors up the generators 
with the scales of the area forms transverse to them). 
In our case, two of these three elements survive: there 
are natural invariant definitions of the generators and of 
the strong conformal geometry, but the conformal struc- 
ture on the space of generators is not preserved under the 
identifications. 

A priori, while for each value of r\ each generator of 
3 + (n) has an affine structure, it is not evident that there 
is a preferred way of identifying these for different values 
of 77. However, because the strong conformal geometry 
and the space of generators are well-defined, the identi- 
fications will extend in a natural way to vectors tangent 
to the generators. This means that supertranslations are 
well-defined. We shall see that there is a natural way 
of measuring the supertranslation relating u m = to 
u V2 = 0, and this will allow us to appropriately account 
for the change of section of the twistor space when we 
compare angular momentum. 

The failure of an invariant conformal structure to exist 
on the space of generators (that is, to be preserved un- 
der evolution) turns out to mean that that in comparing 
the energy-momenta at different cuts higher- j represen- 
tations appear. 

To see this, let us first recall that, since the space of 
generators has naturally the smooth structure of an ori- 
ented sphere, a conformal structure on it is equivalent to 
a complex structure. In the Bondi- Sachs case, the tran- 
sition functions (Cm. ^ (Cm 1 Cm ) preserve this com- 
plex structure, and are thus fractional linear transforma- 
tions. On the other hand, these fractional linear trans- 
formations are isomorphic to the (proper, isochronous) 
Lorentz transformations. Thus in the Bondi-Sachs case, 
the admissible coordinate changes induce Lorentz trans- 
formations on the space of generators. One builds up the 



spaces of asymptotically constant spinors, vectors, etc., 
as functions on the space of generators, and it is the fact 
that the conformal transformations induce Lorentz mo- 
tions which is responsible for these fields breaking neatly 
into Lorentz-invariant representations. 

Now let us turn to the present, non-Bondi-Sachs 
case. If we compute the component P(r]i;( m ,C m ) 01 
the energy-momentum at r\\ along the null vector deter- 
mined by the Bondi stereographic coordinates (C^nCm.) 
in the chart at r\\, we find, as usual, that, as a function 
of (Cm i Cm) 1 the energy-momentum consists of j = 
and j = 1 components, forming a covector. However, 
if we want to compare the energy-momenta at 771 and 
j]2, we must express them both in a common chart, say 
the chart (Cm>, ( V2 ) on the space of generators for C + (772 ) ■ 

Then of course Pi^'Xm^m) wm have only j = and 
j = 1 components, but, because the change of variables 
(C>?i > Cm ) ^ {Cm ' Cm ) wm n °t b e a fractional linear trans- 
formation, the expression for P{f]i) in terms of (C^ 2 , C^ 2 ) 
will generally contain not just j = and j = 1 com- 
ponents, but those for all integral j. So comparison of 
energy-momenta at different surfaces §(r/i), §(%), will 
require higher-j representations. (This sort of behavior 
occurs even strictly at null infinity for angular momen- 
tum, but is a finite-size effect for energy-momentum.) 

Again, this potential failure of a linear identification 
of asymptotic covectors (and elements of the spin-tensor 
algebra generally) as 77 changes is a finite-size effect; it 
will become negligible if one takes the family §(77) of ex- 
traction surfaces distant enough (that is, close enough 
to the physical space-time's null infinity) , assuming that 
the system is indeed isolated. Thus the occurrence of 
these nonlinearities in a numerical computation would be 
a signal that the extraction surface had not been taken 
large enough that a model null infinity, with the usual 
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Bondi-Sachs structure, stable under evolution existed. 

Some further, more technical, discussion of structure 
is given in Subsection IVIIEI 

How the Identifications are Made 

A few words now about how the identifications of J + (77) 
for different 77 are made. They all grow out of two con- 
siderations, which we have already used extensively. The 
first of these is that at any point p(ji) G §(77), there is a 
null geodesic orthogonally outwards in M(if) whose end- 
point lies on 3 + (77) ; holding r\ fixed but varying p{rj) we 
get a preferred cut C(?7)of 3 + (rf). The second is that 
at any p(ji) £ §(77), there is a canonical isomorphism 
Tp(Mphys) — T p (M{rj)) between the tangent space of the 
physical space-time and that of the auxiliary space M(rf). 
(Cf. footnote [17].) 

We define the identifications at the infinitesimal level 
(that is, for infinitesimally separated 77), and then inte- 
grate. At the infinitesimal level, the problem comes down 
to understanding how the the structures on J + (r) + drj) 
are represented by quantities on 3 + (rj). 

Let us consider a one-parameter family of points 77(77) G 
§(77) (with 77 the parameter) in M p h ys - At each point we 
have the null normal l a orthogonally outwards form §(77) , 
and its associated null geodesic in M p h ys - As 77 varies, the 
vector itp hys connecting this family of geodesies is a Ja- 
cobi field, which in turn is determined by its initial data 
(Mp hys , M phys ) at §(77). On the other hand, we may use 
the isomorphism of tangent spaces to regard these as ini- 
tial data for a Jacobi field u a in M(n). We will take the 
limiting value of this field at 3 + (r)) as the definition of 
the rate of change of the end-point of the null geodesic 
orthogonally outwards from p{rj) as 77 varies. The impor- 
tant thing to note here is that this vector, while it lies 
in the tangent space to 3 + (i]), represents the end-point 
of a geodesic with base-point p(r\ + drj) for an infinitesi- 
mally differing value of 77. This is the root of all of the 
identifications. 

We may apply this in several ways. If, for example, we 
require the curve 75(77) to be such that the corresponding 
vector at 3 + (r]) points along a generator, we may say the 
generator of null infinity does not change with 77(77) ; this 
gives the identification of the generators of null infinity 
for different values of 77. Or if we imagine a congruence of 
curves, say p(£, C, 77) with stereographic coordinate we 
get a vector at each point of the cut of 3 + (77) labeled by 
Q, and this vector field over the cut gives a measure the 
supertranslation induced by changing 77 infinitesimally, 
that is, in passing from §(77) to §(77 + drf). 

Outline of this Section 

Subsection IVIIAI gives the computation of the Jacobi 
fields, and Subsection IVH 51 the main formulas for com- 
paring the null infinities. Then Subsections IVII CIVII Dl 



give the formulas for treating the evolution of the energy- 
momentum and angular momentum. The final subsec- 
tion discusses some technical aspects of the structure of 
the null infinities; these are not needed for the computa- 
tions in this paper but are given for completeness of the 
conceptual framework. 

Coordinates 

In what follows, we will be comparing structure on 
the timelike hypersurface T with that on 3 + (rj), and also 
structure on 3 + (?7i) and J + (772) ■ As a ready reference, 
here is a summary of the coordinate systems to be used. 

Recall that, for each 77, we have already defined coor- 
dinates (C,C) = (CvXv) on &(v)- We may regard (C, C) 
then as defined over the whole of T. 

We shall eventually use coordinates (77, (, £) on T. 
However, in the next subsection it will be convenient 
to briefly use coordinates (??,?,£) where £ need not be 
simply related to Q. 

As already indicated, we will have Bondi coordinates 
(C^Cj^r,) on J + (77). 

A. The Jacobi Fields 

We have a family of spacelike surfaces §(77) forming a 
timelike hypersurface T, with v a V a r) > for any future- 
pointing vector v a tangent to T. In practice, it is conve- 
nient to represent the evolution from one surface to the 
next by a connecting vector field w a , that is, a field tan- 
gent to T with w a V a rj — 1. The freedom in choosing w a 
is the freedom to add a vector field which, at each 77, is 
tangent to 8(77), that is, is a linear combination of m a 
and m a . We shall see below that there is a natural way 
to fix this freedom, but for now we leave it unspecified. 

It will be helpful to briefly use coordinates adapted to 
the foliation of T by 77 and the integral curves of w a . 
Near any point of interest on T , fix ?y and let (£, £) be a 
complex coordinatization of §(770). (The use of a complex 
coordinate is only for brevity of treatment; the coordinate 
£ need not be a holomorphic function of £, or have any 
other special relation to it.) We may let the integral 
curve of w a with coordinates (£,£) at 770 be p(r], £,£). If 
we Lie-drag £ along w a (so w a V a S, = 0), then (77, £, £) 
provides a coordinatization of a portion of T. In these 
coordinates, we have w a = d/drj, and d/d£ is a linear 
combination of m a and m a . 



1. Definition of the Fields; their Initial Data 

In this subsubsection, we work in the physical space- 
time, to avoid cumbersome notation, however, the Jacobi 
field is denoted simply u a , rather than Wp^. 

Along each integral curve f>(?7, £, £) of w a (holding £ 
fixed), let 7(77, s,£, £) be the affinely parameterized null 
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geodesic outwards from §(77) in M p h ys , so that l a = 
d s 7{v,s,^,I) and 7(77, 0,£,£) = p(v,£,0- (The scale of 
the afhne parameterization will not matter.) Then differ- 
entiating along 77 we get a Jacobi field u a = d v j(ri, s, £, £) 
connecting these geodesies. 

It is this Jacobi field u a we wish to work out; more 
precisely, we wish to work out u a modulo terms propor- 
tional to l a . The field is determined by its initial data. 



One of those data is simply u a 



= w a ; the other is 



l b V b u a 



One constraint on this second datum is that 



s=0 



lJ b VbU a = 0. (This follows by differentiating l a l a = 0: 
we have = u b X? b {l a l a ) = 2l a u b V b l a = 2lJ b \7 b l a .) The 
other constraint (affecting one complex degree of free- 
dom) comes from requiring that the geodesies meet §(77) 
orthogonally. 

The condition that the geodesies 7(77, s, £, £) meet §(77) 
orthogonally is that (d^p) a l a =0. If we differentiate this 
along w a , and we apply the conditions w b \/ b (d^p) a = 
(d^p) b \7 b w a (which holds since w a = d/dr) in the (77, £, £) 
coordinates on T) and w b V b l a = u b V b l a = l b X7 b u a 
(which holds because u a is a connecting vector for the 
geodesies with tangents l a ), we find 



((d (P ) b v b w a )i a + (d (P )j b v b u a = . 



(71) 



However, since (d^p) b is a complex basis vector span- 
ning the tangent space to §(77), and so the equation is 
equivalent to one with this vector replaced by any other 
complex basis vector for this tangent space; in particular, 
it is equivalent to 



{m b V b w a )l a + m a l b V b u a = , 
which determines the complex datum m a l b, V b u a 



(72) 



s=0 



We shall not need the coordinate £ in what follows. 



2. Solving the Jacobi Equation 

We have derived the initial data for the Jacobi fields 
on the physical space-time M p h ys - We now make use 
of the isomorphism of tangent spaces at §(77) to regard 
these same data as determining Jacobi fields in the aux- 
iliary space-time M(rf), and solve the Jacobi equation 
there. (Since this isomorphism is Wp hys 1— > u a and we 
have already dropped the "phys" subscript, this amounts 
to simply using the formulas derived above for their data 
at 8(77).) Ultimately we are interested in the vectors on 
3 + (77) determined by the asymptotic forms of the Jacobi 
fields. 

In this subsubsection, we work in a frame parallel- 
propagated along the null geodesies, and express the Ja- 
cobi fields in terms of their initial data. In the next one, 
we will transform the asymptotic form of the Jacobi field 
to the Bondi frame. 

It will be convenient to put 



u 00 l a + u 01 



u 10 m" 



u n 'n a 



(73) 



Then the constraints we have worked out above will be 
the initial conditions for the geodesic deviation equation: 



s=0 



u; a and l a V aU 1 



s=0 



l b m a V a w b ■ (74) 



recall that w a is the vector field connecting §(77) to 
§(7/ + drj), so the quantities on the right-hand sides in 
Eq. (174]) are known. As noted above, that the Jacobi 
field represent a null geodesic entails additionally 



l a V a u 11 ' =0. 



(75) 



The geodesic deviation equation itself (l a V a l b V b u c — 
l p l q R prq c u r ) becomes in terms of the components 

7'i 00 ' = fxii 01 ' + * 2 " 11 ' + conjugate (76) 
= -*i77 n ' (77) 
= 0, (78) 



I 10 ' 
•.11' 



where the dots are differentiation with respect to s. (Here 
of course the components of u a and also the curvatures 
^1, ^2 are evaluated at points along the geodesic, not 
at S = §(77); we temporarily violate the convention that 
quantities unsubscripted by NU or B are evaluated at 
S.) Integrating these with the initial conditions (using 
the explicit form of V^i provided by integrating ^ using 
©), we find 



.11' 



00 



,11' 



(79) 



and 
■ W '(s) 



u 



w 10 ' + sl b m a V a w b + 7/; 11 '*i(4|cr| 2 )- 1 x 

{-log(l-a(p+M))(l-a(p-M)-2p« 

- \a\-\p - M)(l - s(p + \a\) log(-*(p + H) 

+ \a\-\p+\a\)(l-s(p-\a\)log(-s(p-\a\)}, 

(80) 



where, on the right-hand side, the spin-coefficients and 
the curvature 'i'l (as well as the field w a ) are evaluated 
at §(77) (we restore the convention about subscripting). 
We note the asymptotic behavior 



,io' 



s(l b m a V a w b + w u 'x) , 



(81) 



where 



X = * x ( -2-V|- a p + 4-V|-V - |<r| 2 )log 



(In the limit \a\ < \p\, we have X w -(3/9) _1 *i.) 



(82) 



We will not need 



,00' 



3. Transformation to the Bondi Frame 

In the previous subsubsection, we found the Jacobi 
field (modulo terms tangent to the geodesic) in a parallel- 
propagated frame. We here transform the asymptotic 
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form of this field to the Bondi frame previously con- 
structed for M(rj). 
We have 

u a modulo l a = u m 'm a + u w 'm a + u 11 ' n a . (83) 



Now m a differs from 



'NU 



by a term proportional to l a , 



and so we may replace m a by m^u in this expression. 
As s — > oo, we may also, to leading order, replace n a 
by ngfy. To see this, first note that we have seen that 
m 01 ',u 10 ' = O(s), u 11 ' = 0(1) ass^oo. Onthe other 
hand, we have n a = — Qm a — Qm a + QQl a . Since 
Q is O(l), making this substitution in Eq. (j8"3")) would 



only change the coefficients of m a , m a (or m^, 
by subdominant terms. Thus 



u a modulo l" 



u 01 'm^u + u w 'fff NU + u n 'nNu 



NU 



-(z^v^ + ^'x)^- 1 ^)^- 

-{l b m a V a w b + w 11 ' X)p- 1 {5'0^ r 



(84) 



Here d/du v = n^, where is a Bondi retarded time 
coordinate for 'J + (rj) adapted to the frame defined by 
(C„,C„) (andP B ). 

Equation (|8"4"f represents the displacement, in 3 + {rf), 
of the cut formed from the null vectors orthogonally out- 
wards, as the two-surface moves from §(77) to §(77 + drf) 
along the vector field w a . It thus codes the relation be- 
tween the null infinities J + (r]) and J + (77 + drf); our next 
task is to develop this into formulas for transition func- 
tions. 



B. Comparison of Null Infinities 

We now have the tools to compare the null infinities 
J + (?7) associated with different values of 77. We recall that 
(for each rj) the invariant structures of J + (77) are its space 
of generators, the conformal structure on that space, and 
the "strong conformal geometry;" we shall see that the 
first and last of these can be identified under changes 
of 77, but not the conformal structure on the space of 
generators. 



1. Identification of the Space of Generators 

We have a family of two-surfaces §(77), and for each 
of these the null geodesies outwards determine a cut of 
the corresponding null infinity We saw in the last 

subsection, however, that we could represent the cut of an 
§(77 + ^77) infinitesimally perturbed from §(77) by a vector 



field in CJ + (77) . Precisely, if w a was a vector field at §(77) 
with w a V a r] = 1, then Eq. ([54]) gave the corresponding 
apparent displacement of the cut. 

We may use this to identify the generators of 'J + (r)) 
for different values of 77. The connecting field w a will 
preserve the generator if the corresponding field at J + (77) 
points purely up the generator, which, from Eq. (|84l) . is 
evidently if 



I'm^aWb+w 11 ' X = 



Expanding l b m a V a wi ) in spin-coefficients we find 



,.11' 



,.01' 



pw 



.10' 



w n 'x. 



(85) 



(86) 



Here w 11 = w a l a depends only on the displacement of 
8(77 + ^77) relative to §(77); it is insensitive to the horizon- 
tal components of w a , which are w 01 and its conjugate. 
We may thus use (1861) as an equation to determine the 
horizontal components of w a from the condition that (|8"Tj) 
vanishes. After a little algebra, we find 



.10 



(p 2 - w\ 



,11' 



,11 



Xw 
-Xw 



IX' 



IX' 



(87) 

We may therefore determine a vector field w a on T by 
requiring w a V a i] — 1 and its components tangential to 
8(77) to be given by Eq. (|57|) . Each integral curve of this 
vector field corresponds to a generator of null infinity, in 
the sense that under the identification of the J + (77) for dif- 
ferent tj's described here, the null geodesies orthogonally 
outwards from §(rj) along this curve are all considered to 
strike the same generator. 

An equivalent way of expressing this is in terms of 
transition functions for the angular coordinates. Let us 
write Cr;o for the stereographic coordinate on J + (t7o) de- 
termined by restricting the stereographic coordinate £ to 
§(?7o) (and identifying the cut of 3 + (r]o) with S(?7o) by 
using the ideal end-points of the null geodesies orthogo- 
nally outwards). We may then extend Cno t° J + ( 7 7) for 
all 77 by requiring 7« a V a C»7 = 

Now let 3 (Co j CotVo, v) be the value of £ = at §(77) 
corresponding to the same generator as does the value Co 
at S(?7o). Then 3 (Co, Co, Vo, v) = Cv C^ 1 can be regarded 
as the transition function from Cn to C> ; , with 77, ?y pa- 
rameterizing the particular choices of coordinate function 
of interest. 

The derivative of 3 with respect to 77 will be the com- 
ponent tip of w a , in the coordinates (77, C>C), if W<1 i s 
chosen to preserve the generators of null infinity. It will 
be conceptually useful to put w a — + tug, where 
Tijy = w 00 l a + w 11 n a are the "vertical" components and 
lug = m a + w w m a are the "horizontal" components 
(with respect to the foliation of T by 77 and the in- 
duced metric). Then = w a V a C — w $ + w h where 
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iu" h V a C- We have 



= w 10 m a + conjugate 

= -w 10 's(pM a +WM a ) + conjugate 

= -w 10 's(pP— +aP-^=) + conjugate 

d 

= —s(w w p + w 01 a)P— + conjugate 

/ / - d 

= -(w 10 p + w al a)p~ 1 S(— + conjugate, 



and so io£ = -(w 10 ' p + w 01 ' a)p^ 1 SC Using Eq. (J57}, 
this becomes = — (-qui 11 + Xw 11 )p~ 1 S(, and so the 
transition function 3 is determined by 

$■ = -ffiw 11 ' +Xw 11 ')p- 1 5( + w$ (89a) 
3(Co,Co,?7o,f7o) = Co- (89b) 

(This has been written as a family of ordinary differential 
equations parameterized by Co, Co, V to emphasize that 
it is a simple evolution equation; one could write d/dr] 
in place of d/ drj just as correctly.) 

In Eqn. (|89aj) . the two terms on the right represent 
contributions which are, respectively, intrinsic and ex- 
trinsic to S (77) in their dependence on the coordinate C- 
The first term can be thought of as representing the rate 
of deformation of the coordinate C due to changes in the 
geometry of the surface as 77 is increased; the second may 
contribute to this but also takes into account the free- 
dom in specifying C on successive surfaces. (Recall that 
C is determined only up to a fractional linear transfor- 
mation by the geometry; it was suggested the remaining 
freedom in C be fixed by comparison with the numerical 
coordinate system in order to have as direct an interpre- 
tation as possible. Cf. the paragraph following the one 
containing Eqn. (|2"3)l .) 



2. The Conformal Structure of the Generators 

As emphasized above, because of finite-size effects, the 
conformal structure of the space of generators is not pre- 
served under evolution. The question of how severe this 
issue is governs the degree to which there is a sense 
of asymptotically constant spinors, vectors and tensors 
which is stable under evolution. This subsubsection dis- 
cusses the obstruction. 

The space of generators has naturally the structure of 
a smooth oriented two-sphere, and a conformal structure 
on this is equivalent to a choice of complex structure, 
and this in turn is equivalent to giving a complex basis 
vector (or covector) up to proportionality. In out case, 
we could take the basis vector to be M a ; then L w M a , 
modulo terms proportional to M a , would give a measure 
of the rate of change in complex structure. Equivalently, 



the shear M a L w M a of M a along w a would measure the 
rate of change. This quantity can be readily computed 
but it is not directly useful here. 

The more direct way of accounting for the change in 
complex structure is the part of d^/drj (Eqn. (|89a)l 1 an- 
tiholomorphic in C; were d^/dr/ holomorphic, the com- 
plex structure would be unchanged at first order in 77 
and di/dr) would induce an infinitesimal Lorentz motion. 

In practice, it is likely that d^/dr/ = will be close 
to holomorphic. We saw above that there are two con- 
tributions to it, one — — (3u> n + Xw 11 )p~ 1 5( de- 
pending on C intrinsic to §(77) and one w§ depending on 
the extension of C off §(rj). For the intrinsic one, we have 
X w — (3p)~ 1 '$>i, and so, if S(i]) is in fact within the 
peeling regime we will have X ~ 0(R~ 3 ). If the sphere 
is nearly round and w 11 is nearly constant, then will 
be small. 

The quantity represents the rate of change of the 
stereographic coordinate as one moves normal to 8(77) 
along w a . This means that w§ depends on just how C 
is extended off §(f]), which in turn depends on how the 
underlying numerical coordinates extend to the future 
of §(rf). As typically the extraction surfaces are very 
distant, nearly round, and these features are reflected to 
good approximation in the numerical coordinates (and 
preserved under evolution), we expect that will give 
something which is close to an infinitesimal fractional 
linear transformation. 



3. The Strong Conformal Geometry 

The "strong conformal geometry" of null infinity links 
the scales of vectors along the generators with the scales 
of those transverse to the generators Q . It can be char- 
acterized by the quantity 



(27)- 1 |P B |- 2 rfCArfC 



d_ 

du ' 



(90) 



where the forms are defined on the space of tangent vec- 
tors to 3 + (r]) modulo the vectors tangent to the genera- 
tors, and the square root is taken in the sense of a linc- 
bundle-valued quantity. (The root taken is irrelevant; 
there are also other equivalent ways of characterizing the 
structure.) This quantity is independent of the Bondi 
frame. 

We have so far constructed, for each 3 + (rj), a stereo- 
graphic coordinate Cij on the space of its generators, and 
we have transition functions C^C^ 1 relating these for dif- 
ferent values of r\. On any 3 + (rj) we may define a Bondi 
coordinate u v with respect to the Bondi frame defined 
by Cjj (and Pb), fixing the zero of u v to lie on the pre- 
ferred cut. Thus the quantity (|90|) is well-defined, and 
requiring it to be preserved under changes of 77 leads to 
a transformation law for the vectors d/du^. 
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We have 



(2 l )- 1 \P B \- 2 dC v A d( v = (2z)- 1 (l + ICr,! 2 )- 1 ^, A d( v 



The full transformation law for the Bondi retarded 
times will be, from Eqn. (f9"2")l again and this, 



_ (l+ICol 2 ) 2 

- (i+ICl*)* 
and so we have 



d 

du Tl 



(2i)- 1 (l+|C J7 o| 2 )- 1 ^oA^ D) (91) 
1 + ICI 2 



+ Au n2 (r ]2 , m ). (94) 



1 + |Cr 12 



■5 r] I 



'/o 



(9 



9u 



(92) 



'»!0 



^. Identification of the Generators 

We now turn to the supertranslations identifying the 
generators of 3 + (r]) for different values of r\. 

We saw above (Eqn. ([S3))) that |A| 2 w fc Z b n° represents 
the apparent displacement of the cut corresponding to 
S(?7 + drf) with respect to the preferred cut in J + (r/). 
To integrate these infinitesimal displacements, however, 
we must express them all as elements of the same vec- 
tor space. We may do this by using the transformation 
rule (f92")) relating n° = d/du n to = d/du m derived 
above; the infinitesimal displacement, expressed as a vec- 



One can verify directly that these transition functions are 
compatible, that is, computing u„ s either directly from 
u m , or from u m in terms of u ril , gives the same answer. 

What we shall actually need is to refer the angular mo- 
menta at different values of 77 to a single value 770- We 



therefore define u ; 

itly, 



u„ o u. 



Ill) 



u(»7i) 



I/O 



1 + ig 2 
1 + iCmi 2 



»?o=0 



Au^(?7, rjo). Explic- 



\X\ 2 w a l a d V . 



(95) 



In this integral, we have ( n = 3(Co, Co, Vo, v) m order 
to keep the generator of null infinity fixed. This ap- 
tor at J+(?7o) for some fixed 770, is (1 + \( Vo | 2 ) -1 (1 + P nes not om y to trie explicit factors of ( v , ( Vl , but 

also to the dependences \X\ 2 w a l a and the derivative 
terms. (So dC, m /dC, n = (d( v / d^)' 1 (dC, v J dC, V0 ) = 

(ds(Co, Co,Vo, v)/d(o) (#3 (Co, Co, Vi, »?)/%)•) 
C. Evolution of Energy Momentum 



\( ri \ 2 )\d( ri0 /d( ri \\X\ 2 w b lbn^ . Thus the supertranslation 
taking the cut labeled by 771 to that labeled by 772 will 
be, in ^+(770), 



Au w (772,?7i) 



1 + lGj 



9C 



I/O 



I^V*^. (93) 



In this integral, the coordinates (Crj ,Cj)o) are ne ^ fixed 
and C„ = 3(CijtnCi7o)- 



To compare the energy-momentum at §(77) with that 
at §(770), then, we express P^AyKA^A 1 by giving tta' as a 
function of C; we then insert for this function 3. We have 



P aa 'w a tta> = P 00 V o7 rcy + P al 'Tf TTv + P W 'wnr > + P u '7fi7iy 

= (i/2)(P°°' + p"') + (i/2)(P°°' - ^ n ')^ T |p - p01 'ytW - pW 'iTW (96) 



In this formula, the components of P a are evaluated at 
§(77) as in Section V. The formula here gives the compo- 
nent of the energy-momentum at §(77) in the direction 
specified by C at S(?7o)- 

As noted above, if the extraction surfaces 8(77) are 
far enough away that they provide good models of cuts 
of null infinity, then 3 will be a fractional linear trans- 
formation representing a Lorentz transformation and 
P Tatta' will be interpretable as the component of 
a covector along the null direction specified by C- An- 
other way of saying this is that energy-momentum (|96[) 
will have only j = and j — 1 components. Because of 
finite-size effects, however, we cannot expect this to hold 
exactly, and there is a question of principle of how to 



extract the (co)vectorial part of the energy-momentum 
when these effects cannot be neglected. 

The natural thing to do is find the boost relative to 
which the energy-momentum (|96[) has zero dipole mo- 
ment (J = 1 component), project the energy-momentum 
in this frame (that is, keep only the j = Q component 
in this frame), and then boost back to the asymptotic 
laboratory frame at §(770). (There will be a unique frame 
in which the dipole moment is zero [15].) While there 
is no simple closed-form expression for this, there is an 
iterative procedure which one would expect to converge 
rapidly. 

The projection of the energy-momentum relative to 
the frame defined by a unit future-pointing vector t AA 
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will be 



E. Two Technical Points 



P AA '(t) =2 5 / 2 7T- 1 J 



ICJ 2 c 
C i 



P BB 'n B ir B , 



(i-MV)l*Cl s 



(97) 



(f0> _£f01' _01O' + |(|2 t ll')3 



-d§. 



where P BB tt b -kb' is given by (|96p . Thus if, starting from 



any t AA , we define 



Ci = OP™ (^P 11 (t n )-\P 01 {t n )\ 2 )- l P AA (i„) (98) 

then the sequence i^" 4 will converge to the time-direction 
determined by the Bondi-Sachs energy-momentum, with 
P AA (t n j converging to that energy-momentum. If 
3(C> C) i s close to £ (as would often be expected), then 
it would be natural to choose = P a / \J PbP b . It may 
well not even be necessary to move to further t® in the 
sequence to attain the accuracy required in many situa- 
tions. 



D. Evolution of Angular Momentum 



The formulas derived above for the evolution of the 
energy-momentum and angular momentum were the 
present paper's goals. As emphasized above, these results 
include possible finite-size corrections, which, if signifi- 
cant, should be interpreted as signs that the extraction 
surfaces are not distant enough to give a stable model of 
null infinity. 

I mention here two further issues related to these finite- 
size effects, points which do not figure in the results above 
but would be relevant if one were to try to draw broader 
lessons for the development of quasilocal kinematics from 
these results. 

The first is that in general the comparisons of energy- 
momentum and angular momentum at S(^o) and §(?7i) 
depend not just on these surfaces themselves but on the 
intermediate ones §(77), t)q < 77 < 771. In other words, 
because of finite-size effects, one would not expect an 
integrable comparison. This issue, of course, would dis- 
appear if the surfaces were actually at null infinity. 

The second issue is that the structures discussed here 
do not actually determine how to evolve the phase of a 
twistor or spinor. This does not lead to any difficulties 
in the formalism given here, but it would be a point to 
keep in mind in developing a more general theory. 



The strategy for comparing the angular momenta at 
different cuts is similar to that for comparing the energy- 
momenta. The essential difference is that we must refer 
all angular momenta to the same "origin" (that is, the 
same cross-section of J + ). We have already found that u 
is the supertranslation relating the measurement of the 
angular momentum at §(77) to that at 8(770). Since a su- 
pertranslation acts on the twistors by simply being added 
to Ab in the parameterization (|67p , (|68p , in order to refer 
the angular momentum back to the original cut 6(770), 
we need to replace Ab by Ab — u in Eq. ([o^)l . as well as 
replacing £ by 3. 

Since, even at a single cut, the angular momentum in 
general relativity is given by an object with components 
for arbitrary j > 1, the angular momentum cannot be re- 
duced to a vectorial object and so the sort of projection 
procedure which was used for the energy-momentum is 
not needed. On the other hand, if it is desired to com- 
pute the components of /i for different values of j, this 
can certainly be done. However, there is no simple al- 
gebraic (that is, involving only finitely many operations 
for each term) transformation taking a resolution of /x in 
spherical harmonics at §(77) (such as might be found via 
Eq. (|7D1) and Tables ED and transporting it to §{r] ), 
even if 3(£, £) is given by a fractional linear transforma- 
tion, because the asymptotic laboratory frames at §(77) 
and §(770) might be relatively boosted, and boosts mix 
infinitely many j- values. (Of course, if 3 can be approxi- 
mated as differing from the identity only to a finite order, 
then an algebraic transformation can be derived.) 



VIII. USERS' GUIDE 

The preceding sections have covered the derivations of 
formulas for the energy-momentum, angular momentum, 
and comparisons of them at different times. The aim of 
the present section is to give a users' guide to the results. 

The starting-point is a spacelike two-surface § of spher- 
ical topology (or, for evolution, a one-parameter family 
§(77) of such surfaces, with V a ?7 timelike). We assume 
that a null tetrad adapted to § has been chosen, and the 
Newman-Penrose quantities at § are available in terms 
of this tetrad. 

The first step is to find a complex stereographic co- 
ordinate C on §; see the paragraph containing Eq. (pM| 
and the two paragraphs thereafter. With this known, 
the factor A giving the rescaling to a Bondi-Sachs frame 
is determined by Eq. (|26[) . This defines an asymptotic 
laboratory frame. 

The second step is to compute the angular potential 
Ab for the Bondi shear. This may be done either via a 
Green's function (using Eqs. <[S5|). ((3"Tj) . (j3"2"10 . or by res- 
olution in spin- weighted spherical harmonics (eqs. Q34[) . 
((33)) , ((ST)]) ; see also the last paragraph of Section Mil for 
phase conventions). 

The third step is to compute the components A a p of 
the kinematic twistor. These are given by integrals over 
S. Table ITvl lists the integrands, which require, besides 
the quantities already discussed, the asymptotic forms 
^1, ^3 °f the Weyl curvature components in the 
poor man's no-incoming-radiation approximation; these 
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asymptotic forms are given in Table U in terms of quan- 
tities on §. (As noted in the text, Table U lists these 
under the assumption \<r\ <C \p\ on §. which should be 
very good for most purposes. Section II shows how to 
compute them more accurately, if required.) 

With the components of A a p known, the Bondi-Sachs 
energy-momentum may be read off directly in the asymp- 
totic laboratory frame: see Eq. . 

The angular momentum is reported as a function 
M(C> C) of the asymptotic direction with respect to asymp- 
totic reference frame, with /_t + ~fi giving the energy- 
moment in the direction and i[i — i]I the spatial angular 
momentum about that axis. One could choose to either- 
present this function directly or to give its resolution in 
spherical harmonics. For a direct presentation, the func- 
tion is given by Eq. (|69p . and this can be resolved into 
spherical harmonics by standard means. If the compo- 
nents of Ab in spherical harmonics have already been 
computed, then the resolution of (i is given by Eq. (|70| . 
which makes use of Tables [V] and I VII 

If the spin and center-of-mass are required, one 
must transform to a boosted asymptotic frame in 
which the time axis lies along the Bondi-Sachs energy- 
momentum [9j. 

To study the evolution of the energy-momentum and 
the angular momentum, one must, besides computing 
them on the different surfaces §(77), also give an invariant 
method for relating the quantities on one surface to those 
on another. If w a is a vector from §(77) to §(77 + drj), one 
solves Eqns. (|5§)) . ([M]) . to find the functions £, rf) and 
u (C;Cj 7 /0 expressing the appearance of the cut 6(8(7/)) 
relative to 6(8(770)). Then Eq. (|96]) (with P AA ' defined 
by (|59D evaluated at §(77)) gives the energy-momentum 
at §(rj) in the directions as specified by the angular vari- 
ables at S(?7o)- This function may, because of finite-size 
effects, in general not be simply a vector but will have 
components in all j > representations. To project the 
vectorial part invariantly, one solves (|9"5)) , (|9"T)) iteratively. 
The angular momentum /1 at §(77) but referred to the 
Bondi coordinates constructed at §(770) is given by re- 
placing Ab by Ab — u and £ by 3 in Eq. There is 
no simple formula for the evolution of the components 
of the angular momentum in spherical harmonics in the 
most general case; one must compute the components 
from the evolved /j, (cf. Section VII D). 



Appendix: Some Properties of Spin- Weighted 
Spherical Harmonics 



This paper relies on some technical properties of spin- 
weighted spherical harmonics, which are derived here. 
The conventions are those of ref. fl4j . 



Definitions 

The spin- weighted spherical harmonics are determined 
as follows. Fix a spin frame o A , i A (normalized with 
6 A i A = 1 and with 2 1 / 2 t AA ' = d A d A ' + i A i A '). Put 

Z{j, m) B ...CD...E = 0(s • ■ • o c Id • ■ • i>E) ■ (A- 1 ) 

j—m j+m 

Now let o A , l a be a second normalized frame, which is 
considered to vary and to determine a point on the sphere 
(corresponding to the null vector o A o A , say). Then one 
puts 



s ^j,m 



Z(j,m) B ...Eoi 



C D 
■ O L 



j+s 



Y — (— 1V+ 171 Z 



3,"i 



(2j + l)!(2i)! 



47r(j + s)!(j - s)!(j + m)!(j - m)! 



(A.2) 



(A.3) 



One has s Y^ m 



(-1) 



— s 1 j,—m • 



The spin-weighted spherical harmonics as defined 
above are functions on certain line bundles over the 
sphere. However, it is common to represent them by 
ordinary functions, by giving their values on preferred 
sections. There are two main conventions for this. In 
the first, the spin- frame is adapted to the complex stere- 
ographic coordinate £ and given by 

o A (CX) = t(l + \(\ 2 y 1/2 (-Cd A -i A ) (A.4) 

l a (CX) = *(i + |C| 2 )- 1/2 M A + C^); (A.5) 

in the second, the adaptation is to the polar coordinates 
6, (f>, and 

o A (9,(f>) = e l ^ /2 cos{0/2)o A + e-^/ 2 sm(9/2)i A (A.6) 
t, A (6,ci>) = -e^ 2 sin(e»/2)6 A + e-^/ 2 cos(0/2)i: A . (A.7) 

These frames differ by a phase only; one has 6 A (6, (f>) = 
-e-^/ 2 6 A (C,C) = -(C/C) 1/2 o A (C,C)- Note that this 
means that, viewed as ordinary functions, we have 



s^j',m (^5 ' 



(-(C/0 1/2 ) 2 %^n(C,C) 



(A. 



In this paper, the spin- weighted harmonics are applied 
at two stages. First, they are used in the solution of the 
twistor equation; in this case, the variable frame is deter- 
mined by oJ| (and the fixed frame by the asymptotic ref- 
erence frame) . The second occurrence of the harmonics is 
in parameterizing the itA' spinors appearing in the def- 
initions of energy-momentum and angular momentum. 
In those case, it is tta which takes on the role of the vari- 
able spinor oa in the spherical harmonics. We note that 
in this case 



7T0' = TTAO = -1/2^1/2,-1/2 
VTl' = TTAL A = --1/2^1/2,1/2 



(A.9) 
(A.10) 
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and similarly 



Tensor Products 



7ro'7r ' 
7T0'7Ti' 

7T0'7fl 
7r '7f ' - -KyTfi 
7Ti'7f 



_iZi,i = v^Ts -i^m 

--1^1,0 = -V 471 "/ 6 -1*1,0 
= v^73-i^i,-i 
0^1,1 = V 47r/6 0^1,1 
2 o-^i,o = 
—qZi i 



(A.11) 
(A.12) 
(A.13) 
(A. 14) 

v / 4V3 n,o (A.15) 
= -V^?6 0^1,-1 ■ (A.16) 



Behavior Under Inversion 

Each of the spin-weighted spherical harmonics is de- 
fined as a function of the spinor oa- We consider here 
how the harmonics change when the spinor is acted on 
by a spatial inversion 

There is a choice of sign in lifting the inversion from 
vectors to spinors; we use ola i-> 2 x l 2 tAA'ol A . (Then the 
spinor tta< appearing in the treatment of angular momen- 
tum is the image of tta 1 under inversion.) 

Each of the harmonics is given as a (normalization fac- 
tor times a) function Za...cd...fo a ■ ■ ■ o c i D ■ ■ ■ l f , where 
there are j + s omicrons and j — s iotas. The antipodal 
map gives oa > — i>a, i— * oa- This will evidently effect 
a change s Yj tTn n> (— iy +s - s Yj, m . (Reversing the sign 
in the definition of the antipodal map on spinors would 
change the action on spin-weighted spherical harmonics 
to s Yj_ m i-> (— l) J_s _ s lj ;TO , that is, would contribute an 
extra minus sign for half-integral spin weights. In this pa- 
per, since all final quantities have integral spin-weights, 
the sign convention, as long as it is kept fixed, is unim- 
portant.) 

However, when one represents the harmonics by ordi- 
nary functions, their behavior under inversions appears 
more complicated, because the sections used to effect the 
trivialization of the bundles are not invariant under in- 
versions. Indeed, evaluating the sections at the antipodal 

point (whose stereographic coordinate is — ( ), we find 

o A (-V,-C l ) = -(C/C) v V(C,0, (A.17) 

^(-r'.-c 1 ) = (c/o 1/2 o A (cx)- (A.i8) 

Thus 

aiw-r'.-c 1 ) - (-i) j+s (c/cr -^(co, 

(A.19) 

where the extra factor is the conversion from the gauge 
at one point to its antipodal point. 

For the frame adapted to the polar system, one has 

o A (ir -6,cf) + n) = -u A (9,<p) (A.20) 
l a (it -6,4> + ir) = -io A (8,(f>), (A.21) 



and so 



We derive here the resolutions of certain products of 
spin- weighted spherical harmonics ("Clebsch-Gordan de- 
compositions") which are used in the text. Because in 
all cases one of the factors has small values of j and s 
(indeed, j = 0, 1, s = —1,0,1) it is easiest to proceed 
iteratively. 

Our starting-point is the identity 



OLaPb...DE = Ol(A0B...DE) - + ■^ € A(B^C...E)FOt t , 

(A.23) 

where cca is any spinor and (3b. ..de is any totally sym- 
metric spinor of valence 2j. Taking a a = oa or a = la 
and /3 B . ..e = Z(j,m) c ...E, we find 



(A.24) 



Z(l/2,-l/2) A Z(j,m) B ...E 

=Z(j + l/2,m- 1/2) ab...e 

7 ~\~ 777 

+ ^p^eA( B Z(j - 1/2, m - l/2) c ... E) 

and 

Z(l/2,l/2) A Z(j,m) B ...E 

=Z{j + l/2,m+ 1/2) ab...e 

~ 'kjT\ tA{BZ{3 ~ 1/2 ' m + x i^ c - e ) ■ 

Contracting now with either o A or t A , and with o B • ■ • o c 
(j + s times) and i D ■ ■ ■ i E (j — s times) , we find 

1/2-^1/2, ±1/2 sZj, m = g+l/2^7+l/2,m±l/2 

tiTm)(J-s) (A.26) 

T TTTTTTrT- , 7~\ s+l/2^j-l/2,m±l/2 



(A.25) 



(2j)(2j + l) 



and 



-1/2^1/2, ±1/2 sZj.m — s-l/2-^j±l/2,m±l/2 

(j T m)(j + s ) (A.27) 
± (2j)(2j + l) -V2^-i/2, m ±i/2 ■ 

We may by repeated application of these build up all 
the tensor decompositions. The cases we need are as 
follows. 

We have -iZ lfi = -1/2^1/2,1/2-1/2^1/2,-1/2- Using 
this, we find, after some algebra 



-\Z\fi \Zj.m —oZj+i lT n — oZj m 



(f-rntW + T) (A - 28) 
1/N.r 1: j - 1,ro ' 

Similarly -\Z\ t ±\ = {-i/2 Z \/2,±i/i) 2 , lrom which 

777 ,J_±rn 

-1^1, ±1 l^j.rn — 0^j+l,m±l ± — — 0^j,m±l 

{j + T m){j T m - I) 

1 ; ~ -o ~ 0^j-l,m±l 



sYj^ji^TT 



. + 7T) = H) 2 ->_ S Y}, 



(A.22) 



4j(4j 2 - 1) 



(A.29) 
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Using oZi,±i = 1/2^1/2^1/2 -1/2^1/2, ±1/2 we find (1/2) -1/2^1/2,1/2 -1/2^1/2,-1/2, 

7 7 _ 7 bT^Q'Tm-l) ,-2_„2 

4(4/^-1) 0^1,0 0-^j,m — O^j+l.m T , 2 _ , 0^j-l,m ■ (A. ,51) 

(A.30) - H J ~ J 

and, using ^i,o = (I/ 2 ) 1/2^1/2,1/2 -1/2^1/2,-1/2 + 
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